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ABSTRACT 

Using data from the mid-infrared to millimeter wavelengths for individual galaxies and for stacked 
ensembles at 0.5 < z < 2, we derive robust estimates of dust masses (Mdust) for main sequence (MS) 
galaxies, which obey a tight correlation between star formation rate (SFR) and stellar mass 
and for star-bursting galaxies that fall outside that relation. Exploiting the correlation of gas to 
dust mass with metallicity (Mgas/Mdust— ■^), we use our measurements to constrain the gas content, 
CO-to-H2 conversion factors (aco) and star formation efficiencies (SEE) of these distant galaxies. 
Using large statistical samples, we confirm that aco and SEE are an order of magnitude higher and 
lower, respectively, in MS galaxies at high redshifts compared to the values of local galaxies with 
equivalently high infrared luminosities (Lir > 10^^ L©). For galaxies within the MS, we show that 
the variations of specific star formation rates (sSFR=SFR/M,) are driven by varying gas fractions. 
For relatively massive galaxies like those in our samples, we show that the hardness of the radiation 
field, {U), which is proportional to the dust mass weighted luminosity (im/Mdust), and the primary 
parameter defining the shape of the IR-spectral energy distribution (SED), is equivalent to SFE/Z. 
For MS galaxies with stellar mass log(Af,/Af0) > 9.7 we measure this quantity, (U), showing that it 
does not depend significantly on either the stellar mass or the sSFR. This is explained as a simple 
consequence of the existing correlations between SFR— M,, — Z and Afgas— SFR. Instead, we show 
that {U) (or equally iiR/Afdust) does evolve, with MS galaxies having harder radiation fields and thus 
warmer temperatures as redshift increases from z = to 2, a trend which can also be understood 
based on the redshift evolution of the A/» — Z and SFR— A^* relations. These results motivate the 
construction of a universal set of SED templates for MS galaxies, that are independent of their sSFR 
or Af,, but which vary as a function of redshift with only one parameter, (U). 

Subject headings: 



1. INTRODUCTION 

Deep and wide multi- wavelength extragalactic surveys 
have greatly enhanced our understanding of galaxy evo- 
lution. It has now been well established that the star for- 
mation rates (SFRs) in galaxies were on average higher 
in the past, with galaxies emitting the bulk of their bolo- 
metric energy in the infrared (IR) progressively domi- 
nating the star formation density of the Universe that 
peaked at z > 1 (e.g., Le Borgne et al. 2009). 

A recent major step forward in characterizing the na- 
ture of star formation in distant galaxies has been the 
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discovery that the majority of star-forming galaxies at 
every redshift define a narrow locus in the stellar mass 
(Af*) — star formation rate (SFR) plane. This correla- 
tion has been observed in the local Universe (Brinchmann 
et al. 2004; Peng et al. 2010), as well as at intermediate 
redshifts 0.5 < z < 3 (Noeske et al. 2007; Elbaz et al. 
2007; Daddi et al. 2007; Pannella et al. 2009; Rodighiero 
et al. 2010a, Karim et al. 2011, Magdis et al. 2010), and 
beyond (e.g., Daddi et al. 2009, Stark et al. 2010), with 
a normalization factor that increases rapidly with look 
back time, mirroring the increase in the star formation 
activity of the early galaxies. 

The fact that this correlation between SFR and 
appears to be present at all epochs has been used to de- 
fine a characteristic specific star formation rate (sSFR = 
SFR/ A/,) at each redshift and stellar mass, and from 
its evolution with cosmic time, a main sequence (MS) 
mode of star formation that is followed by the majority 
of star-forming galaxies. This main sequence also serves 
as a tool to distinguish starburst (SB) from normal^^ 
galaxies at any redshift, simply by measuring the excess 
sSFR of the galaxy with respect to the SFR-A^* corre- 
lation at that redshift. Such outliers are known to exist 
at any redshift. In the local Universe, a class of galaxies 
that are usually identified as starbursts are the (Ultra) 
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Luminous Infrared Galaxies (ULIRGs, Ljr > 10^^ L0, 
Sanders & Mirabel 1996, Elbaz et al. 2007). However, in 
the distant Universe it has proven to be harder than ini- 
tially thought to associate pure starbursts with a specific 
class of high— 2: galaxies, although at least a fraction of 
Submillimeter Galaxies (SMGs) seem to exhibit elevated 
specific star formation rates (Tacconi ct al. 2006; 2008; 
Daddi et al. 2007; 2009, Takagi et al. 2008). 

The remarkable uniformity of MS galaxies, as indicated 
by the SFR-M, correlation, and their contrasting nature 
to that of starburst galaxies, has recently been further 
manifested with respect to two more sets of observable 
parameters. Using direct measurements of the total in- 
frared luminosity (iiR=i8_iooo/im) of distant galaxies 
that have now become possible with the Herschel Space 
Observatory (Pilbratt et al. 2010), Elbaz et al. (2011) 
showed that galaxies that follow the SFR-M* correlation 
are also part of a secondary, infrared main sequence, de- 
fined by a universal total-to-mid-IR luminosity ratio, IRS 
= I/ir/Ls) where is the rest-frame S jim luminosity. 
The majority of MS sequence galaxies, at all luminosi- 
ties and redshifts, arc found to obey this linear correla- 
tion between Lm and Lg, suggesting that they share a 
common IR spectral energy distribution (SED), with a 
mid-to-far-IR shape that has evolved little with cosmic 
time. Interestingly, starburst galaxies, which systemat- 
ically fall above the SFR-M* main sequence, are also 
outliers to the LjR-Lg relation, exhibiting elevated IRS 
values with respect to normal galaxies. 

The third quantity that highlights the distinct nature 
of MS and starbursts encompasses the main driver of star 
formation activity: the molecular gas mass (Mpja), i.e., 
the raw material out of which galaxies form stars. Recent 
studies, by Daddi et al. (2010) and Genzel et al. (2010), 
found that normal galaxies at any redshift have lower star 
formation efficiencies (SFE = SFR/Mhs), compared to 
star-bursting systems that exhibit an accelerated mode 
of star formation activity, probably triggered by a ma- 
jor merger. These studies all suggest that there are two 
regimes of star formation: i) a long-lasting mode, fol- 
lowed by the majority of galaxies at every redshift that 
also form a tight MS in the SFR— and in the Lir — Lg 
planes, and ii) a short -lived starburst mode for galax- 
ies that can become strong outliers from both relations. 
However, this picture heavily relies on estimates 
which are still a matter of significant debate due to the 
poorly determined conversion factor to derive molecu- 
lar gas mass from CO luminosities (aco= My^^/L'qq), 
which is known to vary as a function of metallicity and, 
perhaps, of the intensity of the radiation field (e.g., Leroy 
et al. 2011, Magdis et al. 2011b). 

Although the separation between starbursts and MS 
galaxies is already apparent in the observed Lir- 
L'qq plane, the distinct nature of the star formation ac- 
tivity of the two populations becomes truly evident when 
different aco factors are used for the two classes of galax- 
ies. This differentiation of the aco value has clearly 
been seen in the nearby universe, with local ULIRGs 
(starbursts) having, a value of aco that is a factor of 
~ 6 smaller than that for local spiral galaxies^^ (e.g., 
Downes & Solomon 1998), and seems to be true at high 

We note that Papadopoulos et al. (2012) presented evidence 
that higher aco values are possible in local ULIRGs 



redshift too. Using a dynamical and a dust-to-gas mass 
ratio analysis of star-forming disks at z ^ 1.5, Daddi 
et al. (2010a) and Magdis et al. (2011b) argued for a 
CO conversion factor^^ aco ^ 4.0 for this kind of ob- 
jects, similar to that of the Milky Way and local spi- 
rals. On the other hand, several studies place an upper 
limit of aco ^ O.S for several SMGs (e.g., Tacconi et al. 
200s, Carilh et al. 2010, Magdis et al. 2011b, Hodge et al. 
2012). Furthermore, recent numerical simulations indi- 
cate a clear distinction between the aco values for disks 
and mergers at all redshifts, although with a consider- 
able scatter (e.g., Narayanan et al. 2011, Feldmann et al. 
2012). These findings challenge the common approach 
of "blindly" applying a (local) ULIRG-like aco = O.S 
value to derive the of high— ^ ULIRGs. Instead, 

they highlight the necessity of determining aco in larger 
samples of high— z galaxies and of investigating how it 
varies between normal galaxies and mergers/ starburst, 
but also among galaxies on the main sequence. 

To this end, in Magdis et al. (2011b), we applied, for 
the first time at high redshift, a method for measur- 
ing aco that is commonly used in the local universe. 
The method relies on measuring the total dust mass of 
a galaxy (Mjust) and assuming that it is proportional to 
Mgas (e.g., Leroy et al. 2011). Taking advantage of the 
detailed characterization of the peak of the SED enabled 
by Herschel data, as well as of the Rayleigh- Jeans tail 
from ground based mm observations, we acquired robust 
Mdust estimates for a normal disk at z ~ 1.5 and a star- 
bursting SMG at z = 4.05. The dust mass estimates were 
subsequently used to infer an aco value of ~ 4.0 for the 
disk and an upper limit of ^ 1 for the starburst galaxy, 
in agreement with previous dynamical estimates (Daddi 
et al 2010a, Hodge ct al. 2012). This study also offered a 
first hint of a close link between the SED shape, as traced 
by the dust- ghted luminosity {Lib,/ M^ust) and aco and 
SFE, with MS galaxies exhibiting lower iiR/Mdust val- 
ues compared to SB galaxies. 

Having demonstrated the applicability of the method 
at high redshift, here we extend this exercise to a larger 
sample of main sequence disks at z ~ 1.5 and z ~ 0.5 
that have extensive photometry from rest-frame UV to 
mm wavelengths, including state of the art Herschel data 
from the Herschel Great Observatories Origins Deep Sur- 
vey {GOODS- Herschel, PI D. Elbaz) as well as low- 
transition CO observations and mm continuum interfero- 
mctric measurements. We use this to investigate the av- 
erage aco value for high— z MS galaxies, finding consis- 
tent values with what is predictcxi l)y the local aco ver- 
sus Z trend and what is inferred by dynamical analyses 
(Daddi et al. 2010b). 

Pushing the method further, we here analyse large sta- 
tistical samples of galaxies at z ^ 1 and z ^ 2 to derive 
their average dust and gas mass content. While the lack 
of CO measurements for such large samples prevent us 
from learning an}rthing on their conversion factor, the 
recovered information is invaluable to study the distri- 
bution of gas in high redshift galaxies. This allows us to 
investigate possible variations of gas fractions and SFE 
within the MS, providing insights about the origin of the 

The units of aco, Mqpc^'^ (K km s^^)^^, are omitted from 
the text for brevity. Also, aco estimates in this work account for 
the presence of Helium coexisting with the molecular hydrogen 
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"thickness" of the SFR-M* relation that has been shown 
to reflect the variation of the physical properties of MS 
galaxies (e.g., Salmi et al. 2012, Elbaz et al. 2011). Fur- 
thermore, the available data allow us to investigate the 
SED shape of MS galaxies as function of redshift and 
also as a function of their their offset from the main se- 
quence. Put together, in this study we attempt to charac- 
terize the gas, dust and SED properties of main sequence 
galaxies throughout cosmic time and therefore provide a 
better understanding of the nature of the sources that 
dominate the star formation density at all redshifts (e.g., 
Rodighiero et al. 2011, Sargent et al. 2012). 

This paper is organized as follows. In Section 2, we 
describe our sample and the multi-wavelength observa- 
tions. In Section 3 we present a detailed analysis of the 
far-IR properties. In Section 4 we discuss the method 
to derive aco and present aco a-nd SEE estimates for 
MS and starburst galaxies. In Section 5, we explore the 
variations of SFE within the MS, and discuss two pos- 
sible scenarios to explain the observed dispersion of the 
SFR-M* correlation. In Section 6 we investigate these 
two scenarios throiigh stacking, while in Section 7 we ex- 
plore the shape of the SED of MS galaxies as a function 
of cosmic time and build template SEDs of MS galaxies 
in various redshift bins. Finally, in Section 8, we pro- 
vide a discussion motivated by the results of this study, 
while a summary of the latter is presented in Section 9. 
Throughout the paper we adopt fim = 0.3, Hq = 71 km 
s^i Mpc~\ Oa = 0.7 and a Chabrier IMF (Chabrier et 
al. 2003). 

2. OBSERVATIONS AND SAMPLE SELECTION 

The aim of this study is to investigate the far-IR and 
gas properties of high— z, normal, main sequence galax- 
ies. Instead of using a large, heterogeneous sample, we 
decided to focus on a small but well defined set of sources 
selected to meet the following criteria: 1) available spec- 
troscopic reclshifts; 2) rich rest-frame UV to mid-IR pho- 
tometry; 3) no excess in the specific star formation rate, 
i.e., sSFR/sSFRms < 3 (where sSFRms is the MS trend, 
as detailed below) 4) available [1-0] or [2-1] low-transition 
CO observations, to enable CO luminosity estimates, 
without the caveat of the uncertainties introduced by 
the excitation corrections affecting higher CO transition 
lines; 5) Herschel detections that offer a detailed sam- 
pling of the far-IR part of the SED; and, if possible, 6) 
mm continuum data that provide a proper characteriza- 
tion of the Rayleigh- Jeans tail, which is crucial for robust 
dust mass estimates. 

To derive a characteristic sSFRms at a given redshift 
and a given stellar mass, we define a main sequence, 
SFRms(-2) M^,), that varies with stellar mass with a slope 
of 0.81 (e.g., Rodrighero et al. 2011), and evolves with 
time as (l + z)^-^^ (e.g., Elbaz et al. 2011, PanncUa et al. 
2009). In what follows we describe the Herschel observa- 
tions used in this study, as well as the sample of galaxies 
considered here. 

2.1. Herschel Data 

We use deep 100 and 160 PACS and 250, 350, and 
500 Atm SPIRE observations from the GOODS- ffersc/ie/ 

program. Details about the observations are given in El- 
baz et al. (2011). Herschel fluxes are derived from point- 
spread function (PSF) fitting using galfit (Peng et al. 



2002). A very extensive set of priors, including all galax- 
ies detected in the ultra-deep Spitzer Multiband Imaging 
Photometer (MIPS) 24 /xm imaging, is used for source ex- 
traction and photometry at 100, 160 and 250 /um, which 
effectively allow us to obtain robust flux estimates for 
relatively isolated sources, even beyond formal confusion 
limits at 250 /im. For 350 and 500 /im, this approach does 
not allow accurate measurements due to the increasingly 
large PSFs. Hence, we use a reduced set of priors based 
primarily on Very Large Array (VLA) radio detections, 
resulting in flux uncertainties consistent with the confu- 
sion noise at the IR wavelengths. Our measurements are 
in good agreement with the alternative catalogs used in 
Elbaz et al. (2011). The advantage of using galfit for 
PSF fitting is in its detailed treatment of the covariance 
matrix to estimate error bars in the flux measurements, 
which is crucial to eventually derive reliable estimate of 
flux errors for the case of blended/ neighbouring sources. 
The effective flux errors in each band can vary substan- 
tially with position, depending on the local density of 
prior sources over areas comparable to the PSF. A de- 
tailed description of the flux measurements and Monte 
Carlo (MC) derivations of the uncertainties will be pre- 
sented elsewhere (E. Daddi et al., in preparation). We 
correct the PACS photometry for a small flux bias intro- 
duced (due to source filtering for background subtrac- 
tion) during data reduction (see H-GOODS public data 
rclease^^; Popesso et al. in preparation). 

2.2. A Sample of Main Sequence Galaxies at z 1.5 
and z ~ 0.5 

Daddi et al. (2010) presented PdBI C0[2-l] emission 
line detections of six star- forming galaxies at z ~ 1.5, 
originally selected by using the "star-forming BzK" color 
criterion (Daddi et al. 2004b). All six sources (BzK- 
21000, BzK-17999, BzK-4171, BzK-16000, BzK-12591 
and BzK-22536) have spectroscopic redshifts (now con- 
firmed by multiple CO detections; for optical redshifts 
see Stern et al. in preparation, Cowie et al. 2004) and 
robust PACS and/or SPIRE detections. Four of these 
sources also benefit from 1.3 mm continuum data derived 
as a by-product from the CO [5-4] emission line observa- 
tions (Dannerbauer et al. 2012 in prep). For BzK-21000, 
continuum detections and uper limits are also obtained 
at 1.1, 2.2, and 3.3 mm (Daddi et al. 2010a; Danner- 
bauer et al. 2009). In addition to IRAC and MIPS 24 /um 
data, the sources are seen in the 16 /im InfraRcd Spec- 
trograph peak-up image (Teplitz et al. 2011). Although 
we will revisit (and confirm) the star formation rates of 
these sources based on the Herschel data, the existing 
UV, mid-IR and radio SFR estimates, along with stel- 
lar mass measurements derived by fitting the Bruzual & 
Chariot (2003) model SEDs to their rest-frame UV to 
near-IR photometry (Daddi et 2010a), consistently place 
them in the SFR— M* main sequence at 2: ~ 1.5 (Daddi 
et al. 2010). Finally, the UV rest-frame morphologies, 
the double-peaked CO profiles, the large spatial extent 
of their CO reservoirs, and the low gas excitation of the 
sources, all provide strong evidence that they are large, 
clumpy, rotating disks (Daddi et al. 2010a). 

We also consider three z ~ 0.5 normal disks with Sersic 
index n < 1.5 and spectroscopic redshifts, for which 

http://hedam.oamp.fr/GOODS-HerscheI/index.php 
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TABLE 1 

PACS, SPIRE AND 1.3mm photometry for MS GALAXIES IN THIS STUDY 



Source 


RA^ 


DEC 


■5100 


S'leo 


S250 


S350 


•S'soo 


5*1. 3mm 




J2000 


J2000 


mJy 


mJy 


mJy 


mJy 


mJy 


mJy 


ID-8049 


188.9751587 


62.1787071 


10.7±0.7 


20.1±1.5 


23.0±1.5 


23.7±5.2 


12.5±4.5 




ID-5819 


189.1774139 


62.1594429 


12.2±0.7 


19.5±1.3 


13.6±1.8 


5.7±6.5 


3.2±7.1 




ID-7691 


188.9502106 


62.1763153 


14.3±0.7 


19.0±1.3 


13.6±1.5 


7.6±6.5 


4.6±7.2 




BzK-4171 


189.1106110 


62.1431656 


3.4±0.4 


10.3±1.1 


15.1±3.1 


7.5±5.6 


7.0±5.8 


0.04±0.14 


BzK-12591 


189.4224091 


62.2142181 


10.2± 0.7 


19.3±1.4 


25.6±2.0 


19.7±4.4 


9.2±5.6 




BzK-25536 


189.3679962 


62.3152504 


0.8±0.7 


3.0±0.9 


5.2±1.2 


0.2±6.1 


5.6±5.6 




BzK-21000 


189.2942505 


62.3762665 


9.1 ±0.5 


17.0±1.4 


25.3±4.4 


20.1±4.7 


11.6±7.5 


0.83±0.36 


BzK-17999 


189.4658661 


62.2556038 


4.9±0.5 


12.5±1.1 


15.3±1.3 


11.1±5.4 


5.9±4.8 


0.30±0.10 


BzK- 16000 


189.1253967 


62.2410736 


1.8±0.5 


4.4±0.7 


10.2±4.4 


9.5±5.2 


3.3±5.9 


0.53±0.13 



Notes; 

1; coordinates are from VLA 1.4 GHz continuum emission (Morrison et al. 2010). The VLA 1.4 GHz map has a resolution of 1.8" 
the typical position accuracy for our sources is 0.20" 



and 



we have C0[2-l] emission line detections (Daddi ct al. 
2010b). Similar to the z ~ 1.5 sample, the sources are 
detected in the PACS and/or SPIRE bands and are part 
of the SFR— main sequence based on their 24 /um- 
derived IR luminosities and SFRs, which are known to 
be robust in this redshift range (Elbaz et al. 2010, 2011). 
The Herschel and millimetre photometry of the main se- 
quence galaxies considered here are summarised in Table 
1. 

2.3. A Sam,ple of High- z SMGs 

As a comparison sample to our MS galaxies, we also 
include in our analysis a small sample of SMGs: GN20, 
SMMJ2135-0102 and HERMES J105751. 1+573027. The 
selection of the targets was driven by the need for 
available multi-wavelength photometry, including Her- 
schel and mm continuum observations, as well as CO [1-0] 
emission line detections. 

GN20 was already studied in Magdis et al. (2011). It 
is one of the best-studied SMGs to date, the most lumi- 
nous and also one of the most distant (z = 4.05, Pope et 
al. 2006, Daddi et al. 2009) in the GOODS-N field. Her- 
schel photometry and the far-IR/mm properties of the 
source have already been presented in detail in Magdis 
et al. (2011b). In brief, the source is detected in all Her- 
schel bands (apart from 100 /im) and in the AzTEC 1.1 
mm map (Perera et al. 2008) while continuum emission 
is also measured at 2.2, 3.3, and 6.6 mm (Carilli et al. 
2011) as well as at 1.4 GHz with the VLA (Morrison et 
al. 2010). Furthermore, Carilli et al. (2010) reported the 
detection of the CO[1-0] and C0[2-l] hues with the VLA, 
and CO [6-5] and CO [5-4] lines with the Plateau de Bure 
Interferometer (PdBI) and the Combined Array for Re- 
search in Millimeter Astronomy (CARMA), respectively. 
A compilation of the photometric data is given in Table 
1 of Magdis et al. (2011b). 

SMMJ2135-0102 (SMM-J2135 hereafter) is a highly 
magnified SMG at z = 2.325, with an amplification 
factor oi jjL = 32.4 ± 4.5, serendipitously discovered by 
Swinbank et al. (2010) behind the cluster MACS J2135- 
01. The source has exquisite photometric coverage from 
rest-frame UV to radio including SPIRE broadband data 
(Swinbank et al. 2010 Table 1 and Ivison et al. 2010 Ta- 
ble 1). Swinbank et al. (2010) also report the detection of 
CO [1-0] and CO [3-2] emission lines using the Zpectrome- 
ter on the Green Bank Telescope and PdBI observations, 
while Ivison et al. (2010), using the SPIRE Fourier Trans- 
form Spectrometer (FTS), presented the detection of the 



[CII]158 /xm cooling line. 

Finally, HERMES J105751. 1+573027 (HSLW-01 here- 
after), is a 77ersc/ieZ/ SPIRE- selected galaxy at z=2.957, 
multiply-lensed (magnification factor ijl = 10.9 + 0.7) by 
a foreground group of galaxies. The source was discov- 
ered in Science Demonstration Phase /fersc/ie//SPIRE 
observations of the Lockman-SWIRE field as part of the 
Herschel Multi-tiered Extragalactic Survey (HerMES; S. 
Oliver et al. 2012). The optical to mm photometry of the 
source is presented in Table 1 of Conley at al. (2011), 
while Ricchcrs ct al. (2011) reports the detection of C0[5- 
4], CO [3-2], and CO [1-0] emission using the PdBI and 
CARMA and the Green Bank Telescope. 

We note that due to the high magnification factor for 
SMM-J2135 and HSLW-01, these sources might not be 
representative of the bulk population of SMGs, tradi- 
tionally selected as 6*850 > 5 mJy. Also, a large number 
of non-lcnscd SMGs with CO [3-2] observations is avail- 
able in the literature. However, we choose not to con- 
sider them in this study, as the uncertain gas excitation 
results in dubious CO [1-0] estimates (e.g., Ivison et al. 
2010, Riechers et al. 2011), which are essential for inves- 
tigating the gas properties of the sources. 



2.4. A Statistical Sample of z ^ 
Galaxies 



I and z-^ 2 MS 



In addition to our samples of individually-detected MS 
galaxies at z ^ 1.5 and z ~ 0.5, we also perform a stack- 
ing analysis to derive the average SED of MS galaxies 
at 2; ~ 1 and 2; ~ 2. Specifically, we use the 2; ~ 1 
GOODS samples of Salmi ct al. (2012; see also Daddi et 
al. 2007 and Pannclla et al. in prep). In order to remove 
bulge-dominated galaxies with low sSFR, we discard ob- 
jects with a Sersic index n > 1.5 (Salmi et al. 2012), as 
measured from the GOODS ACS images (Giavalisco et 
al. 2004) , and consider only 24/xm detected galaxies from 
the GOODS ultra-deep Spzizer imaging. To remove star- 
bursts we proceed as follows: i) from and redshift, we 
compute sSFRmsj i-e., the fiducial mean sSFR expected 
for MS galaxies with a given mass and redshift, derived 
as described at the start of §2 ii) from Herschel data we 
derive Ljr and subsequently SFR estimates, and exclude 
galaxies with measured sSFR/sSFRms > 3. Since not 
all starburst systems are expected to be detected by Her- 
schel we also omit sources with sSFRs/sSFRms > 3 using 
the UV-based, corrected for extinction, SFR estimates by 
Daddi et al. (2007) and the Spitzer estimates of Salmi et 
al. (2012). The 2; ~ 2 sources are drawn from the BzK- 
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selected sample of Daddi ct al. (2007; now extended to 
deeper K-band depths by Pannella ct al. in prep.). Since 
this color-selection technique naturally excludes quies- 
cent galaxies, we only exclude starbursts identified using 
the same method as for the 2; ~ 1 sample. The total 
number of 2: ~ 1 and 2 ~ 2 galaxies considered in the 
stacking analysis is 1569 and 3618, respectively. The av- 
erage SEDs of these two samples are then measured by 
stacking, using various standard techniques depending 
on the wavelength: 

• 16 fim: We use the Teplitz et al. (2011) maps. 
The stacking is performed using the IAS stacking 
library (Bethermin et al. 2010a) and PSF-fitting 

photometry. 

• 24 and 70 /xm: We use the 24 fim and 70 /xm (Prayer 
et al. 2006b) maps of GOODS, the IAS stacking 
library and aperture photometry with parameters 
similar to that in Bethermin et al. (2010a). 

• 100 and 160 /xm: We use the GOODS- ffersc/ie/ 

PACS data (Elbaz et al. 2011), the IAS stacking 
library and PSF-fitting. We apply the appropriate 
flux correction for faint, non masked, sources to the 
PACS stacks (Popesso et al. in preparation). 

• 250, 350 and 500 /xm: We use the GOODS- 
Herschel SPIRE (Elbaz et al. 2011) data and the 

mean of pixels centred on sources as in Bethermin 
et al. (2012a). The bias due to clustering is esti- 
mated to be about 20% at 500 iim and thus neg- 
ligible compared with the statistical imccrtainties. 
This bias is smaller at shorter wavelengths where 
the beam is narrower (Bethermin et al. 2012a). 

• 870 /im and 1100/im: We used the WciB ct al. 
(2009) LABOCA map of eCDFS and the Perera et 
al. (2008) AzTEC map of GOODS-N. Contrary to 
SPIRE, (sub-)mm data are noise-limited, it is thus 
optimal to beam-smooth the map before stacking 
(Vicira ct al. in prep.). We thus apply the Mars- 
den et al. (2009) method to perform our stacking, 
which takes the mean of pixels centered on sources 
in the beam-smoothed map. 

At all wavelengths, we used a bootstrap technique to 
estimate the uncertainties (Jauzac et al. 2011), and the 

mean fluxes measured in the two flelds arc combined 
quadratically to produce a mean SED at z ~ 1 and 2; ~ 2. 
For the 2 ~ 2 normal galaxies we also construct average 
SEDs in three stellar mass bins in the range of 9.7 < 
log(MH./Mo) < 11.2 with a bin size of A log(M*/M0) = 
0.5, as well as in four bins of sSFR/sSFRms over the 
range -0.4 < log(sSFR/sSFRMs) < 0.4, with a bin size 
of Alog(sSFR/sSFRMs) = 0.2. 

3. DERIVATION OF FAR-IR PROPERTIES 

Several key physical properties of distant galaxies, such 
as infrared luminosities {Lin), dust temperatures (Td) 
and dust masses (A'/dust), can be estimated by fltting 
their mid-to-far-IR SEDs with various models and tem- 
plates calibrated in the local Universe. However, the lack 
of sufficient data for a proper characterization of the SED 
of distant galaxies has often limited this kind of analy- 
sis to models suffering from (necessarily) over-simplified 



assumptions and large generalizations. The Spitzer, Her- 
schel, and millimeter data available for the galaxies in our 
high— 2 sample provide thorough photometric sampling 
of their SEDs, allowing the use of more realistic models 
that previously have been applied mainly in the analysis 
of nearby galaxies. Here we consider both the physically- 
motivated Draine & Li (2007) (DL07 hereafter) models, 
as well as the more simplistic, but widely used, modified 
black body model (MBB). 

3.1. The Draine & Li 2007 Model 

We employ the dust models of DL07, which consti- 
tute an update of those developed by Weingartner & 
Draine (2001) and Li & Draine (2001), and which were 
successfully applied to the integrated photometry of the 
Spitzer Nearby Galaxy Survey (SINGS) galaxies (Draine 
et al. 2007). These models describe the interstellar dust 
as a mixture of carbonaceous and amorphous silicate 
grains, whose size distributions are chosen to mimic 
the observed extinction law in the Milky Way (MW), 
the Large Magellanic Cloud (LMC), and the Small 
Magellanic Cloud (SMC) bar region. The properties of 
these grains are parametrized by the PAH index, ^pah, 
defined as the fraction of the dust mass in the form of 
PAH grains. The majority of the dust is supposed to be 
located in the diffuse ISM, heated by a radiation field 
with a constant intensity Umin- A smaller fraction 7 of 
the dust is exposed to starlight with intensities ranging 
from Umin to Umax, representing the dust enclosed in 
photo-dissociation regions (PDRs). Although this PDR 
component contains only a small fraction of the total 
dust mass, in some galaxies it contributes a substantial 
fraction of the total power radiated by the dust. Then, 
according to DL07, the amount of dust dM^ust exposed 
to radiation intensities between U and U + dU can be ex- 
pressed as a combination of a ^-function and a power law: 



dM, 



dust 



dU 



(1 - 7)Mdust^(U - Umin) + 7Mdust 



1 



(1) 

with {Umin < Umax, a 1); -^dust the total dust mass, 7 
the fraction of the dust mass that is associated with the 
power-law part of the starlight intensity distribution, and 
Umin, Umax, characterizing the distribution of starlight 
intensities in the high-intensity regions. 

Following DL07, the spectrum of a galaxy can be de- 
scribed by a linear combination of one stellar component 
approximated by a blackbody with color temperature T* 
= 5000K, and two dust components, one arising from 
dust in the diffuse ISM, heated by a minimum radiation 
field Umin ("diffuse ISM" component), and one from dust 
heated by a power-law distribution of starlight, associ- 
ated with the intense photodissociation regions ( "PDR" 
component). Then, the model emission spectrum of a 
galaxy at distance D is: 



max 



U" 



f-°'i<^' = fi,B,(T,) 



M, 



dust 



47rD2 



(l-7)pi.°'+7P. , (2) 



where fJ* is the solid angle subtended by stellar pho- 
tospheres, = p^J'\qpAK,U^in) and = Pu{qpAH 
■ t^min, f^inax7 «) are the emitted power per unit frequency 
per unit dust mass for dust heated by a single starlight 
intensity f/min, and dust heated by a power-law distri- 
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Fig. 1. — Observed SEDs of ^ ~ 1.5 BzK galaxies, overlaid with the best fit Draine & Li (2007) models (black) and the best fit single- 
temperature modified blackbody (blue). The black dashed line is the DL07 model without the stellar component that is depicted with a 
green dotted-dashed line. Orange dashed and cyan dot-dashed lines show separate contributions of starlight and emission from dust heated 
hy U = Umin (diffuse ISM component) and dust heated by Umin < U < Umax ("PDR" component) respectively. The fitted parameters 
from the best fit Draine & Li (2007) model fits are listed within each panel. 




ID 100 1000 



Fig. 2. — Same as Figure 1, but for the z ~ 0.5 galaxies in our sample. 



bution of starlight intensities dM/dU oc [/ " extending 

from [/„iin to t/max- 

In principle, the dust models in their most gen- 
eral form are dictated by seven free parameters, 
(r2,, qPAH, t^min, C^max, «, 7 and Md). However, Draine et 
al. (2007) showed that the overall fit is insensitive to the 
adopted dust model (MW, LMC and SMC) and the pre- 
cise values of a and Umax ■ In fact they showed that fixed 
values of a = 2 and Umax = 10^ successfully described 
the SEDs of galaxies with a wide range of properties. 
They also favor the choice of MW dust models for which 
a set of models with qpAH ranging from 0.4% to 4.6% 



is available. Furthermore, since small Umin values corre- 
spond to dust temperatures below ~ 15 K that cannot 
be constrained by far-IR photometry alone, in the ab- 
sence of rest-frame submm data, they suggest using 0.7 
< Umin < 25. While this lower cutoff for Umin prevents 
the fit from converging to erroneously large amounts of 
cold dust heated by weak starlight {Umin < 0.7), the price 
to pay is a possible underestimate of the total dust mass if 
large amounts of cold dust are indeed present. However, 
Draine et al. (2007) concluded that omitting submm data 
from the fit increases the scatter of the derived masses 
up to 50% but does not introduce a systematic bias in 
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Fig. 3. — Same as Figure 1, but for the SMGs in our sample. 



the derived total dust masses. 

Under these assumptions, we fit the mid-IR to mm 
data points of each galaxy in our sample, searching for 
the best fit model by minimization and parametrizing 
the goodness of fit by the value of the reduced x^, xt = 

/df (where df is the number of degrees of freedom). 
The best fit model yields a total dust mass {M^ust), Umim 
7 and gpAH while to derive Lir^ estimates we integrate 
the emerging SEDs from 8- to 1000/im: 



L 



IR 



lOOO^im 



Lv{\) X — dA. 



(3) 



Si-L'fn 



A by-product of the best fit model is also the dust 
weighted mean starlight intensity scale factor, (U), de- 
fined as: 



L 



dust 



{U) = 



Po^Idust 

(1 - i)u„ 



1 

min 



for a 



.(4) 

where Pq is the power absorbed per unit dust mass in a 
radiation field with U = 1. Note that essentially {U) is 
proportional to Li^J M^^st, and as we will discuss later 
for the definition of Lir adopted here, i.e. Ls-iooo/jm, 
our data suggest that Pq «125. Uncertainties in Lir and 
Mdust are quantified using Monte Carlo simulations. To 
summarise, for each galaxy a Gaussian random num- 
ber generator was used to create 1000 artificial fiux sets 
from the original fluxes and measurement errors. These 
new data sets were then fitted in the same way, and the 
standard deviation in the new parameters was taken to 
represent the uncertainty in the parameters found from 
the real data set. Best fit values along with their corre- 
sponding uncertainties are listed in Table 2 for all sources 
in our sample. To check for possible contamination of 
the submm continuum broad band photometry from C+ 
(158 /im) emission (e.g.. Small et al. 2011), we repeated 
the fit excluding the affected bands (i.e., 350 /im and 
500 /im for z and z ^ 1.5 — 2.0 respectively without 
noticing any effect in the derived parameters. The best 

^ idust quoted below is similar to Lir, but integrated from to 



fit models along with the observed SEDs of the z ^ 1.5, 
z ^ 0.5 galaxies and SMGs are shown in Figures [U [5] and 
[3] respectively. 

3.2. The Importance of Millimeter Data 

While Herschel data accurately probe the peak of the 
SED in the far-IR emission of distant galaxies, rest-frame 
submm observations (Arest ^ 250 /^m) are necessary to 
sample the Rayleigh- Jeans tail. Since the available far- 
IR photometry for five of our sources is restricted to Her- 
schel observations, it is important to explore possible bi- 
ases or systematics introduced by the lack of mm contin- 
uum data. Given the ever increasing number of distant 
galaxies with Herschel photometry, this investigation will 
also provide guidance for similar studies in the future. To 
assess the significance of adding mm photometry in the 
derivation of the far-IR properties of the galaxies, and 
particularly their dust mass, we repeat the fitting proce- 
dure, this time excluding any data at wavelengths longer 
than rest-frame 200 /im, for the four z ~ 1.5 BzK and the 
three SMGs with available mm continuum observations. 

A comparison between the derived dust masses with 
and without the use of mm data is shown in Figure 
m (left). We find no evidence for a strong systematic 
bias, with an average ratio between the two estimates 
(M^^/M^"™'") = 0.80 ± 0.2. Interestingly, the sources 
with the largest discrepancies are the SMGs, although 
the sample is too small to clearly demonstrate the exis- 
tence of a systematic effect. However, the addition of mm 
data has a noticeable impact on the uncertainties of the 
derived Mdust estimates, which are reduced, on average, 
by a factor of ^ 2. Studies in the local universe reach 
similar conclusions, with the addition that in the absence 
of rest-frame sub-mm data, dust masses tend to be un- 
derestimated for metal-poor galaxies. This could serve 
as an indirect indication that our sample mainly consists 
of metal-rich sources, something that we will also argue 
in Section 4. 

The derived 7 and ?7min are also in broad agreement be- 
tween the two cases. However, we notice a weak system- 
atic bias towards higher {U) values, i.e., stronger mean 
radiation fields, when mm data are considered in the fit, 
reflecting the fact that ([/) oc ijR /A/dust- We conclude 
that such trends in {U) and Afdust suggest that mm data 
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TABLE 2 

Physical properties of individually detected sources, based on Draine & Li 2007 Models 



Source 


^spcc 


Xi^ 


log LiR 


logA/d (DL07) 




7 


qpAH 




Td °- 


/3 " 








Lq 


Mq 




% 


% 




K 




ID-8049 


0.507 


2.1 


11.22±0.02 


8.81±0.13 


2.0 


0.4 


2.50 


2.0 


22±2 


1.5 " 


ID-5819 


0.530 


1.78 


11.26±0.01 


8.22±0.10 


7.0 


2.4 


2.50 


8.8 


31±1 


1.5 " 


ID-7691 


0.637 


0.65 


11.54±0.03 


8.28±0.10 


8.0 


7.6 


2.50 


14.5 


31±2 


1.5 ^ 


BzK-4171 


1.465 


1.56 


11.98±0.04 


8.70±0.09 


15.0 


0.3 


2.50 


15.4 


37±2 


1.5 ^ 


BzK-12591 


1.600 


1.04 


12.44±0.02 


9.09±0.09 


10.0 


7.2 


2.50 


17.6 


37±1 


1.5 


BzK-25536 


1.459 


1.94 


11.46±0.06 


8.52±0.26 


7.0 


0.0 


3.19 


7.0 


33±3 


1.5 " 


BzK-21000 


1.523 


1.12 


12.32±0.01 


9.07±0.06 


10.0 


3.5 


2.50 


13.7 


35±2 


1.4±0.2 


BzK-16000 


1.522 


1.71 


11.87±0.03 


9.11±0.07 


3.0 


5.0 


3.19 


4.7 


30±1 


1.5±0.2 


BzK-17999 


1.414 


0.86 


12.06±0.02 


8.78±0.07 


12.0 


2.7 


2.50 


15.3 


33±1 


1.9±0.2 


GN20 


4.055 


1.44 


13.25±0.06 


9.67±0.06 


25.0 


1.8 


1.20 


29.6 


33±2 


2.1±0.2 


SMM-J2135 " 


2.325 


3.41 


12.30±0.06 


8.85±0.06 


20.0 


1.5 


0.47 


21.9 


30±1 


2.0±0.2 


HSLW-01 " 


2.957 


6.61 


13.13±0.06 


8.93±0.08 


25.0 


30.0 


0.47 


134.0 


39±2 


2.1±0.2 


Stack-zl <i 


0.98 


3.91 


11.21±0.09 


8.14±0.13 


8.0 


2.0 


3.90 


9.7 


32±2 


1.5'' 


Stack-z2 


1.97 


2.01 


11.55±0.06 


8.27±0.11 


12.0 


2.5 


3.19 


15.1 


38±2 


1.5'' 


Notes: 



a: Derived based on single temperature MBB 
b: /3 fixed to this value in the MBB fit 
c; Values corrected for magnification 
d: Best fit values accounting for the redshift distribution of the sample 



can place better constraints on the diffuse ISM emis- 
sion, as well as on the relative contribution of a PDR 
component to the total radiation field, and consequently 
on the Mfiust of individual high— z galaxies. Finally, we 
note that the detailed sampling of the peak of the far-IR 
emission provided by Herschel data can result in robust 
Lib. estimates without the need of mm data. 

3.3. Comparison With Modified Black Body Fits 

Another method for deriving estimates for the dust 
masses and other far-IR properties such as dust temper- 
atures and dust emissivity indices is to fit the far-IR 
to submm SED of the galaxies with a single-temperature 
modified blackbody, expressed as: 

U « (5) 

e '=^<i - 1 

where Td is the effective dust temperature and /? is the 
effective dust emissivity index. Then, from the best fit 
model, one can estimate Af^ust from the relation: 

with Krost = i^aiir^)^ 



(1 + z)nrcstBu{Xrest-,Td) ' 



(6) 

where Si, is the observed flux density, Di^ is the lumi- 
nosity distance, and Krest is the rest-frame dust mass 
absorption coefficient at the observed wavelength. While 
this is a rather simplistic approach, mainly adopted due 
to the lack of sufficient sampling of the SED of distant 
galaxies, it has been one of the most widely-used meth- 
ods in the literature. Therefore, an analysis based on 
MBB-models not only provides estimates of the effective 
dust temperature and dust emissivity of the galaxies in 
our sample, but also a valuable comparison between dust 
masses inferred with the MBB and DL07 methods. 

We fit the standard form of a modified black body, 
considering observed data points with \rest > 60 /im, 
to avoid emission from very small grains that dominate 
at shorter wavelengths. For the cases where mm data 
are available we let /3 vary as a free parameter, while 
for the rest we have assumed a fixed value of /3 = 1.5, 
typical of disk-like, main sequence galaxies (Magdis et al. 



2011b, Elbaz et al. 2011). From the best fit model, we 
then estimate the total Mdust through equation 6. For 
consistency with the DL07 models we adopt a value of 
K250tim= 5.1 cm^ g~^ (Li & Draine 2001). To obtain the 
best fit models and the corresponding uncertainties of the 
parameters, we followed the same procedure as for the 
DL07 models. The derived parameters are summarized 
in Table 2 and the best fit models are shown in Figures 
1, 2 and 3. 

A comparison between dust masses derived by DL07 
and MBB models in shown in Figure |4] (right). We see 
that a modified black body appears to infer dust masses 
that are lower than those derived based on DL07 mod- 
els on average by a factor of - 2 {{M^^'^^yM^^^^) = 
1.96 ± 0.5). A similar trend was recently reported by 
Dale et al. (2012), who found that the discrepancy be- 
tween the two dust mass estimates is smaller for sources 
with warmer S'yo/S'ieo colors. Interestingly, when we con- 
volve the best fit DL07 rest-frame SEDs of our galaxies 
with the 70- and 160 fim PACS fihers, we find that the 
sources with the warmest S^q/Siqq are indeed those with 
the best agreement in the dust masses derived from the 
two methods. The reason behind this discrepancy has 
been addressed by Dunne et al. (2000). The grains of 
a given size and material are exposed to different inten- 
sities of the interstellar radiation field and thus attain 
different equilibrium temperatures which will contribute 
differently to the SED. The single-temperature models 
cannot account for this range of Td in the ISM and at- 
tempt to simultaneously fit both the Wien side of the 
grey body, which is dominated by warm dust, as well as 
the Rayleigh-Jean tail, sensitive to colder dust emission. 
The net effect is that the temperatures are driven to- 
wards higher values, consequently resulting in lower dust 
masses. While restricting the fit to A < 100 /im does not 
provide an adequate solution to the problem (Dale et al. 
2012), a two-temperature black body fit is more realis- 
tic and returns dust masses that are larger by a factor of 

2 compared to those derived based on a single Td MBB 
(Dunne & Eales 2001), in line with those inferred by the 
DL07 technique. 

Despite the large uncertainties associated with the 
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Fig. 4. — Left: The significance of the addition of mm continuum data for the derivation of dust masses. A comparison between the 
DL07 M(j„gt estimates inferred with and without the addition of mm data to the fit. Although, the two estimates seem to be overall in 
good agreement, we note that in absence of mm continuum data we tend to overestimate the Afj^gt in SMGs. Right: Comparison between 
dust masses derived based on a single temperature modified black body (MBB) and those derive by the Draine &; Li (2007) models. Filled 
squares and open circles indicate sources with and without mm data. For the former /3 was treated as a free parameter in the MBB fit 
while for the latter it was fixed to /3 = 1.5. Orange stars represent the SMGs considered in this study. In both panels the black solid line 
corresponds to unity and the two grey dashed lines its offset by a factor of 2 and 0.5. 



MBB analysis, it is worth noting that for all SMGs in 
our sample we measure an effective dust emissivity in- 
dex /3 ~ 2.0, in agreement with the findings of Magnelli 
et al. (2012). On the other hand, the far-IR to mm SEDs 
of all BzK galaxies are best described by /3 ~ 1.5, simi- 
lar to that found by Elbaz et al. (2011) for MS galaxies. 
However, we stress that the effective /3 does not neces- 
sarily reflect the intrinsic /3 of the dust grains due to 
the degeneracy between j3 and the temperature distribu- 
tion of dust grains in the ISM. The /3-Td degeneracy also 
prevents a meaningful comparison between the derived 
Td for the galaxies in our sample. 

4. THE GAS CONTENT OF MAIN SEQUENCE AND 
STAR-BURSTING GALAXIES 

Several studies in the local Universe have revealed that 
the total gas-to-dust mass ratio ((5gdr) is correlated with 
the gas-phase oxygen abundance, in the sense that more 
metal-rich galaxies tend to exhibit lower (5gdr values 
(e.g., Muhoz-Mateos et al. 2009, Leroy et al. 2011 and 
references within). In addition to crucial information 
regarding the amount of metals trapped in the dust, 
this correlation can serve as a valuable tool for deriv- 
ing indirect estimates of the CO to molecular gas mass 
(MHa) conversion factor. In particular, if we know the 
metallicity, the dust mass and the mass of the atomic 
gas mass (Mhi) of a galaxy, then we can estimate the 
(5gdr through the Jqdr— relation, and Mjja through 
the following relation: 

<5GDRMdust = Mga. = Mh. + Afei (7) 
Subsequently, if L'^q is measured then we can estimate 
aco , since — aco x . This method has success- 
fully been applied in the local Universe (e.g., Leroy et al. 
2011) and has recently been extended to high— z galax- 
ies by the pilot study of Magdis et al. (2011b). Even 
though restricted to only two galaxies, which are also 
part of this study (GN20 and BzK21000), Magdis et al. 
(2011b) showed that this approach provides aco esti- 
mates that are consistent with other independent mea- 



surements and dynamical arguments. Although it has 
been a common practice to adopt a value of aco = 0-8, 
typical of local ULIRGs, for high— z ULIRGs, our earlier 
study suggested that such low values apply only to the 
subset of genuinely star-bursting systems at high red- 
shift, while high— z main sequence galaxies, even those 
with > 10-^^Lq (i.e., those classified as ULIRGs) 
have larger aco values, similar to those of local spirals. 
Here, we wish to extend this experiment, this time using 
a substantially larger sample, and to investigate possi- 
ble correlations between aco ^^nd starburst versus main 
sequence indicators. 

4.1. Estimating Metallicities 

A key ingredient of this method is the metallicity of the 
galaxies, for which we have to rely on indirect measure- 
ments. We first derive stellar mass estimates by fitting 
the Bruzual & Chariot (2003) model SEDs to their rest- 
frame UV to near-IR spectrum. Then for the z ^ 1.5 
BzK galaxies and SMMJ2135-0102 we use the - Z 
relation at z ~ 2 from Erb et al. (2006), while for the 
z 0.5 sources we adopt the fundamental metallicity 
relation (FMR) of Mannucci et al. (2010) that relates 
the SFR and the stellar mass to metallicity. To check 
for possible systematics, we also derive metallicity es- 
timates for the z 1.5 sources with the FMR rela- 
tion and find that the two methods provide very sim- 
ilar results. For GN20 and HLSW-1 the situation is 
more complicated, as the FMR formula is only appli- 
cable up to z ~ 2.5, and the Erb et al. (2006) relation 
saturates above ^ IO^^Mq. For these two sources 
we adopt the line of reasoning of Magdis et al. (2011b). 
Namely, in addition to the metallicity estimates based 
on the FMR relation, we also consider the extreme case 
where the huge SFR of the two sources 1500 M0 yr~^) 
is due to a final burst of star formation triggered by a 
major merger that will eventually transform the galaxy 
into a massive elliptical. Once star formation ceases, 
the mass and metallicity of the resulting galaxy will 
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Fig. 5. — Left: Mgas/Afdust vs. metallicity for a sample of galaxies in the Local Group by Leroy et al. (2011) (grey circles) and local 
ULIRGs by Downes &; Solomon (1998) (grey stars). The solid black line is the best linear regression fit to Leroy's sample and the grey 
shadowed area depicts the dispersion of the correlation. Filled black squares and orange stars indicate the position of 2 ~ 0.5 — 1.5 MS 
galaxies and high— 2 SMGs respectively, based on the aco vs. metallicity relation derived in the right panel (see eq. 8). The vertical, 
dashed line indicates the limiting metallicity for which the presented relation is used in this study. Right: Inferred «co values, based 
on the 5gdr method, against metallicity for high— 2 MS galaxies (black squares), SMGs (orange stars). Similar to the left panel, we 
include the sample of galaxies in the Local Group by Leroy et al. (2011) (grey circles), as well as the local ULIRGs by Downes & Solomon 
(1998) (grey stars). All metallicities have been calculated on PP04 scale. The blue solid line represents the best regression fit to the data 
(excluding local ULIRGs), with a slope of -1.39. 



not change further, and one might therefore apply the 
mass-metalHcity relation of present-day elliptical galax- 
ies (e.g., Calura et al. 2009). Combining the metallic- 
ity estimates based on this scenario with those derived 
based on the FMR relation, we estimate that the metal- 
hcities could fall in the ranges 12 -h log [O/i?] = 8.8 to 9.2 
for GN20, and 8.6 to 9.0 for HLSW-1. For the analysis 
we will adopt 12 + \og[0/H] = 9.0 ± 0.2 for GN20 and 
12 + log[0/iJ] ^ 8.8 ± 0.2 for HLSW-1. We also adopt 
a typical uncertainty of 0.2 for the rest of the sources. 
The assumed metallicities, all calibrated at the Pettini & 
Pagel (2004) (PP04) scale, along with the stellar masses, 
are given in Table 3. 

4.2. Derivation of the CO to Mfjj Conversion Factor 

To derive aco estimates for our sample, we first derive 
a Squr—Z relation using the local sample presented by 
Leroy et al. (2011), after converting all metallicities to 
the PP04 scale (Figure [5] left). The data yield a tight 
correlation between the two quantities'^ described as: 
log <5gdr= (10.54±1.0) - (0.99±0.12) x {12+\og{0/H)), 
with a scatter of 0.15 dex. Having obtained the 
-^^dust and metallicity estimates, we then use the 
Sgbr—Z relation to derive Mgas and subsequently 
estimate aco from the equation Mgas = c^cox L'^q- 
For the last step we have assumed that at high— z, 
^ Mhi or equivalently that Mgas ~ Mh^, as 
supported by both observational evidence (e.g., Daddi 
et al. 2010; Tacconi et al. 2010, Geach et al. 2011) and 
theoretical arguments (Blitz & Rosolowsky 2006; Bigiel 
et al. 2008, Obreschkow et al. 2009). All values and their 
corresponding uncertainties, which take into account 
both the dispersion of the ^gdr— ^ relation and the 
uncertainties in Z and Mjust, are summarized in Table 
3. In Figure [S] (right), we plot the derived aco values 
as a function of metallicity, along with the estimates of 

Adopting the relation quoted by Leroy ct al. (2011), i.e., 
log (5gdr= (9.4 ±1.1) - (0.85 ±0.13) x (12 + log[0///]), has virtu- 
ally no impact on the results presented here, with a mean difference 
between the derived aco estimates of a factor of ~ 1.15 



Leroy et al. (2011) for a sample of local normal galaxies 
as well as for a sample of local ULIRGs from Downes 
& Solomon (1998), for which we were able to compute 
their metallicities on the PP04 scale. Fitting high— z 
MS galaxies along with the local sample of Leroy et al. 
(2011), yields: 

log(aco) = -(1.39±0.3)x(12+log[0/i/])ppo4+12.8±2.2 

with aco decreasing for higher metallicities, in agree- 
ment with previous studies (e.g., Wilson 1995, Israel et 
al. 1997, Schruba et al. 2012, and references within). It 
is evident that main sequence galaxies at any redshift 
have higher aco values than those of local star-bursting 
ULIRGs and more similar to those of local normal galax- 
ies. In particular, for high— z MS galaxies, we find a 
mean aco of 5.5±0.4, very close to the MW value of 
4.6 (e.g.. Strong & Mattox 1996, Dame et al. 2001). In- 
terestingly, SMM-J2135 falls close to the locus of main 
sequence galaxies with aco ~ 3, while for the remaining 
SMGs we find aco ~ 1, close to the average value of local 
ULIRGs (aco ^ 0.8, e.g., Solomon et al. 1997, Tacconi 
et al. 2008). Finally, simply as a consistency check, we 
derive and plot in Figure [5] (left) the Afgas/Mdust val- 
ues of our targets, using aco estimates as obtained from 
equation 8. 

4.3. Limitations 

Before trying to interpret our findings regarding the 
Q^co values, it is important to discuss possible caveats 
and limitations of the method. A key ingredient for the 
derivation of aco is the Mdust estimate. These estimates 
rely heavily on the adopted model of the grain-size distri- 
bution. Dust masses based on different grain size compo- 
sitions and opacities in the literature can vary even by a 
factor of ~ 3, indicating that the absolute values should 
be treated with caution (e.g., Galliano et al. 2011). How- 
ever, the relative values of dust masses derived based on 
the same assumed dust model should be correct and pro- 
vide a meaningful comparison, as long as the dust has 
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Fig. 6. — iJ^oZ-'^dust against metallicity for the same sample of 
galaxies as in Figure |S] This plot serves as a sanity check for the 
assumption that the A/gas/Afjust vs. metallicity relation does not 
evolve considerably with cosmic time. 

similar properties in all galaxies. Therefore, any trends 
arising from the derived Mdust estimates are essentially 
insensitive to the assumed dust model. 

One important assumption of our technique is that 
the Sgbr—Z relation defined by local galaxies does not 
evolve substantially with cosmic time. These quantities 
are very intimately related, as dust is ultimately made 
of metals. For example, solar metallicity of Zq = 0.02 
means that 2% of the gas is in metals. The local ^gdr— ^ 
relation suggests an average trend where 1% of the gas 
into dust. This means that about 50% of the metals 
are locked into dust. This is already quite an efhcient 
process of metal condensations into dust, and it is hard 
to think how this could be even more efficient at high- 
z (where, e.g., less time is available for the evolution of 
low-mass stars into the AGB phase) , with a much higher 
proportion of metals locked into the dust phase, a situa- 
tion that would lead to an overestimate of Mgaswith our 
technique. Similarly, chemical evolution studies suggest 
a small evolution of dust-to-metal ratio, i.e., by a maxi- 
mum factor of ~ 2 from 2 Gyr after the Big Bang up to 
the present time (e.g., Edmunds 2001, Inoue 2003, Calura 
et al. 2008). Furthermore, if we plot the ratio of the di- 
rect observables, L^Q/Mdust, as a function of metalhcity, 
both local and high— z galaxies follow the same trend, 
suggesting that the local Sqbr—Z relation is valid at 
high— 2; (Figure|6|). Finally, while the Leroy et al. (2011) 
sample yields a slope of —1.0 in the Sgdr—Z relation, 
Muhoz-Mateos et al (2009), find a rather steeper slope 
of -2.45. However, adopting the Mufioz-Mateos relation 
(after rescaling it to PP04 metallicity) , returns very sim- 
ilar aco values, as the two relations intercept at metal- 
licities very close to the average metallicity of our sample 
(12 + log[0/H] w 8.6). We note though that for GN20 
and HLSW (the sources with the highest Z estimates), 
the Mufioz-Mateos relation would reduce aco values by 
a factor of ^ 3. 

4.4. aco o,s a Function of Specific Star Formation Rate 

With robust Lm measurements in hand, we are in a po- 
sition to derive accurate estimates of the SFR and sSFR 
for galaxies in our sample. We first convert the derived 



LiR to star formation rates by using the Kennicutt (1998) 
relation, scaled for a Chabrier IMF for consistency with 
our stellar mass estimates. We then infer the specific star 
formation rate of each source and compare it to the char- 
acteristic sSFRms of the main sequence galaxies at the 
corresponding redshift of the source, using the method- 
ology described in §2 up to z = 2.5. For the two SMGs 
at z ~ 3 and z ~ 4, we use the SFR-M* relations of 
Magdis et al. (2010a) and Daddi et al. (2009) respec- 
tively. In order to conservatively identify MS galaxies, 
we classify starbursts as galaxies with an excess sSFR 
relative to that of the MS galaxies by at least a factor of 
3, i.e., with sSFR/sSFRms > 3. We note however that 
the sSFR/sSFRms indicator is only a statistical measure 
of the star formation mode of the galaxies, and there is no 
rigid limit that separates main sequence from starburst 
galaxies. 

Our analysis suggests that high— z MS galaxies, asso- 
ciated with a secular star formation mode, exhibit higher 
aco values at any redshift compared to those of merger- 
driven local ULIRGs. Since the latter are known to be 
strong outliers from the local SFR-M, relation, we ex- 
pect a dependence of the aco value with sSFR. However, 
the functional dependence between the two quantities is 
not straightforward. In the local universe, there is ev- 
idence for a bimodality in the aco value, linked with 
the two known star formation modes: normal disks (for 
Z ^ Zq) are associated with aco ~ 4.4 and merger 
driven ULIRGs with aco ^ 0.8. As we will discuss later, 
if the aco of normal galaxies depends only on metallic- 
ity, then we expect only a weak increase of aco as a 
function of sSFR/sSFRms within the MS. Consequently, 
in order to incorporate the observed steep decrease of 
aco in starburst systems, aco and the relative offset 
from the main sequence (i.e., sSFR/sSFRms) will be re- 
lated with a step-function. On the other hand, if aco 
also depends on other parameters than metallicity, such 
as the compactness and/or the dumpiness of the ISM of 
normal galaxies, and if the ISM conditions of MS galaxies 
do vary significantly, then aco could smoothly decrease 
as we depart from the MS, and move towards the star- 
burst regime (e.g., Narayanan et al. 2012). 

The scenario of a continuous variation of aco with 
sSFR/sSFRms is depicted in Figure [iK- In addition to 
our galaxies we also include a sample of local normal 
galaxies for which we have robust and SFR estimates 
from Leroy et al. (2008) and aco measurements from 
Schruba et al. (2012). A linear regression fit yields the 
following relation: 

aco = (5.8 ± 2.0) x [sSFR/sSFRms]*""'^^"^"-^^^ (9) 

with a scatter of 0.3 dex. We stress that, given the lack of 
a sufficient number of starbursts in our sample, this rela- 
tion is poorly constrained in the starburst regime. Never- 
theless, the locus of local ULIRGs, seems to be consistent 
with the emerging trend. Note that since we lack a sam- 
ple of local ULIRGs for which both aco and sSFR are 
accurately determined, the position of local ULIRGs in 
Figure [7] indicates the average aco and sSFR of the pop- 
ulation, as derived by Downes & Solomon (1998) and Da 
Cunlia et al. (2010b), respectively. To further check the 
robustness of the fit we also perturbed the original the 
data within the errors and repeated the fit for 1000 reali- 
sations. The mean slope derived with this method is very 
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Fig. 7. — a) Inferred aco values against the offset from the main sequence, a) Inferred aco values for the MS galaxies (black squares) 
and SMGs (orange stars) considered in this study against their offset from the main sequence, parametrized as sSFR/sSFRms- Empty grey 
circles represent local disks from Leroy et al. (2008) and Schruba et al. (2012), while the empty grey star the locus of local ULIRGs. The 
data (only high— z galaxies) are best fit with a slope of -0.67 (red solid line) with a scatter of 0.3 dex and a Spearman's rank correlation 
coefficient of p = —0.51. b) Same as the left panel but for the case of a metallicity-dependent aco- Tracks show the variation predicted 
by Sargent et al. (2012b, in prep) based on (i) the relative importance of the starburst and main-sequence mode of star formation as a 
function of sSFR/sSFRms a^nd (ii) a "metallicity- dependent" acOi where the metallicity of main-sequence galaxies is inferred based on the 
stellar mass— star formation— metallicity relation of Mannucci et al. (2010). Blue (purple), solid, dashed and long-dashed lines correspond 
to z = 1.5 {z = 0.5) for M* = 1.0 X 10^° Mq, 5.0 X 10^° M©, and 1.0 x 10^^ Mq. The grey long-dashed line corresponds to z = 0.0 and 
M, = 1.0 x 10" M0. 



close to the one describing the original data. A Spear- 
man's rank correlation test yields a correlation coefficient 
of p = —0.51, with ap— value of 0.03, suggesting a moder- 
ately significant correlation between sSFR/sSFRms and 
aco- However, repeating the Spearman's test, this time 
excluding star-bursting systems, does not yield a sta- 
tistically significant correlation, with p = —0.15 and a 
p— value of 0.3, indicating a very small dependence, if 
any, between aco and sSFR within the MS. 

It is therefore unclear whether aco varies inside the 
MS sequence, or if instead there is a bifurcation between 
the mean values appropriate for MS versus SB galaxies. 
If the aco conversion factor primarily depends on metal- 
licity, one would expect little or no variations with sSFR 
at constant mass within the MS. Detailed computations 
for this scenario are presented in Sargent et al. (2012b, 
in prep.), in which burst-like activity is assumed to be 
characterised by a constant, low conversion factor, cho- 
sen here to be = 0.8. To predict the variations of 
aco as a function of sSFR/sSFRms according to this 
scenario, Sargent et al. (2012) compute the relative im- 
portance of main-sequence and starburst activity at a 
given position within the stellar mass vs. SFR plane us- 
ing the decomposition of the sSFR-distribution at fixed 
stellar mass derived in Sargent et al. (2012a; based on 
the data of Rodighiero et al. 2011). Metallicities, which 
form the basis for computing aco for an ISM experi- 
encing "normal" star-formation activity, vary smoothly 
as a function of stellar mass and star formation rate ac- 
cording to the calibration of the fundamental metallicity 
plane given in Mannucci et al. (2010; see also Lara-Lopez 
et al. 2010). To compute the corresponding value of the 
conversion factor, Sargent et al. (2012b in prep.) assume 
a variation of aco oc as they ffird that this relation 
between metallicity and aco reproduces the faint-end 
slope of the z = CO-LF of Keres et al. (2003) best. 



under the set of assumptions just described. 

The colored tracks in Figure [TId delineate the varia- 
tions of aco with sSFR/(sSFR)ms expected for galax- 
ies in three bins of total stellar mass. Three regions 
can be clearly distinguished: (i) the main-sequence lo- 
cus, characterised by a gradual increase of aco with 
sSFR that is caused by the slight rise in metallicity pre- 
dicted by the FMR; (ii) the starburst region at high 
sSFR/ (sSFR)ms ^ 4, where the conversion factor, by 
construction, assumes a mass- and redshift-independent, 
constant value; and (iii) a narrow transition region (span- 
ning roughly sSFR/(sSFR)ms £ [2,4]) where aco drops 
from an approximately Milky Way-like conversion factor 
to the starburst value. Note that, in contrast to the pre- 
vious case, aco varies only slightly among MS galaxies. 
This weak dependence of aco on sSFR/sSFRms intro- 
duces a step as we move from MS to starburst systems. 
In the emerging picture, the majority of the star-forming 
population is dominated by either the main-sequence or 
starburst mode, and only a small fraction consists of 
"composite" star-forming galaxies hosting both normal 
and starburst activity. For composite sources, the aco 
should be interpreted as a mass-weighted, average con- 
version factor that reflects the relative amount of the 
molecular gas reservoir that fuels star-forming sites ex- 
periencing burst-like and secular star-formation events, 
respectively. Note also that variations with redshift (in- 
dicated by different colored lines in Figure 7b) are small 
for the limited range of stellar mass covered by our sam- 
ple, and would likely be indistinguishable within the nat- 
ural scatter about the median trends plotted in the Fig- 
ure. 

Both scenarios appear to be consistent with the data, 
leaving open the question of whether the transition from 
normal to starburst galaxies is followed by a smooth or 
a step- like variation of aco- However, both scenarios 
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agree in that aco should not vary much within the main 
sequence. Indeed, even for the case of continuous vari- 
ation, the decrease of aco becomes statistically signifi- 
cant only for strongly star-bursting systems, as equation 
9 indicates that aco would vary at most by a factor 
of ~ 2 within the full range of the MS, well within the 
observed scatter of the relation. Therefore, the aver- 
age value aco ~ 5, as derived in this study, can be re- 
garded as representative for the whole population of MS 
galaxies at any redshift, with stellar masses > lO^^M©. 
We remind the reader that because sSFRms increases 
with redshift as (l-|-z)^-^^, at least out to z w 2.5, galax- 
ies with (U)LIRG-like luminosities can enter the MS at 
higher redshifts or at large stellar masses (e.g., Sargent 
et al. 2012). We stress that our analysis suggests that 
the aco of these systems, i.e., high— z main sequence 
ULIRGs, is on average a factor of ~ 5 larger than that 
of local ULIRGs. 

4.5. Comparison with Theoretical Predictions 

The determination of the aco value has also been 
the focus of several theoretical studies. In particular 
Narayanan et al. (2011), investigated the dependence of 
aco on the galactic environment in numerical simula- 
tions of disk galaxies and galaxy mergers, and reported 
a relationship between aco and the CO surface bright- 
ness of a galaxy. Here, we compare our observationally 
constrained aco values with those predicted by their the- 
oretical approach. According to Narayanan et al. (2011) 



-0.32 



aco 



6.75 X 10^°(Wco) 
6.3 X 1019 X (Z/Zo)o < 



(10) 



where (Wco), is the CO surface brightness of the galaxy 
in K Km s~^. According to Solomon & Vanden Bout 
(2005) 



CO 



23.5 X CIb*s Dl Ico 



where i'co is measured in K km s ^ pc 



(l + z)-3 (11) 
^, ilB*s is the solid 

angle of the source convolved with the telescope beam 
measured in arcsec^, and Iqq is the observed integrated 
line intensity in K km s^^. In the case of an infinitely 
good resolution, as assumed in Narayanan's simulations, 
this becomes the solid angle subtended by the source. 
The CO line luminosity can also be expressed for a source 
of any size in terms of the total line flux: 



Leo = 3.25 X 10^ Sco Av ,.-^1 (1 + z\ 



-3 



(12) 



where 5*00 Av is the velocity integrated flux in Jy km 
s^i. From equations 11 and 12, and for the case of an 
infinitely good resolution, we derive : 



Wco = 1-38 X 10^ 



4 X ly 



obs 



X ScoAv 



(13) 



where fobs is the observed frequency and is the size 
of the source in arcsec^^. For our sources we use, where 
possible, the CO sizes, and UV sizes when CO observa- 
tions are too noisy to reliably measure the source extent. 

We introduce a (1 -|- z) term here, not present in the original 
Narayanan et al. formula, which is required in order to make aco 
independent on redshift for fixed physical properties of the galaxies. 
Without that term the systematic difference between predicted and 
measured aco values would grow to a factor of 2 
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Fig. 8. — Comparison between the <5gdr inferred aco values for 
the galaxies in this study and those predicted by the theoretical 
approach of Narayanan et al. (2011). MS galaxies at 2; ~ 0.5 and 
and 2 ~ 1.5 are shown with black squares and black squares with 
blue circles respectively, while SMGs are depicted as orange stars. 
The solid line is the 1-to-l relation and the dashed line its offset in 
y, by a factor of 2. 

Using equations 10 and 13, we then get an estimate of 
aco for each source in our sample according to the pre- 
scription of Narayanan et al. (2011), and in Figure [5] we 
compare them to the values derived based on our method. 
Although there seems to be a small systematic offset (on 
average by a factor of ~ 1.4) towards lower values for 
the theoretical approach, the overall agreement of the 
aco values is very good when one considers the large 
uncertainties and assumptions of the two independent 
methods. We stress that the derived (VFco) and there- 
fore the estimated aco values are very sensitive to the 
choice of the size indicator. 

4.6. Implications for the Star Formation Activity of 
High—z MS Galaxies 

While the sSFR probes the star formation mode of 
a galaxy only indirectly and from a statistical point of 
view, accurate measurements of its star formation effi- 
ciency, SFE — SFR/Mgas, can provide more direct and 
reliable answers. With the robust aco estimates de- 
rived in this study, we can convert the CO measurements 
into Mgas and subsequently, infer the star formation effi- 
ciency of galaxies in our sample. In what follows we use 
SFE = SFR/Mgas and SFE = Lm/Mgas interchange- 
ably, since SFR is linearly related to Lir through SFR — 
Lin xlO~i° (i.e., through the Kennicutt 1998 relation, 
calibrated for a Chabrier IMF). 

We have so far established that even though there ap- 
pears to be only small variation of aco among MS galax- 
ies, the aco of high—z MS galaxies is similar to that of 
local disks with lower Lir, and is approximately ~ 5 
times larger than the average aco value observed for lo- 
cal ULIRGs. The direct implication of this finding is 
that for a given L'qq, one could expect a large variation 
in the amount of molecular gas for galaxies between MS 
and starburst galaxies. This also suggests that for a given 
Lm. (or equally, SFR) , galaxies in the main sequence have 
lower star formation efficiencies, as compared to star- 
bursting systems. Indeed, using the inferred aco values 
we derive SFE estimates for each source in our sample 
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Fig. 9. — a) Star formation efficiency, SFE=LiR/Afgas , against the offset from tfie main sequence. Black squares and orange stars are 
the MS galaxies and the SMGs considered in this study. Empty grey circles are local normal galaxies (Leroy et al. 2008 and Schruba et al. 
20f2) and grey stars depict local ULtRGs (Zaurin et al. 2009 and Solomon et al. f997). For local ULlRGs we have assumed (aco=0-8)- 
The red solid line depicts the best fit to the data with a slope of 1.34, while a Spearman's rank test indicates a tight correlation between 
the plotted quantities with p = 0.71. b) Same as the left panel but for the case of a weak increase of SFE with sSF'R/sSFRjvrs within the 
main sequence. Tracks show the dependence of SFE on sSFR/sSFRms a-s expected for a slightly sub-linear relation between Afgas and SFR 
and the relative importance of the starburst and main-sequence mode of star formation as a function of sSFR/sSFRms- Blue (purple), 
solid, dashed and long-dashed lines correspond to z = 1.5 (z = 0.5) for A/, = 1.0 X lO"' M©, 5.0 X 10^° M0, and 1.0 X 10^^ Mq. Grey 
long-dashed line corresponds to z = and Af« = 1.0 X 10^^ Mq. 

and plot them against sSFR/sSFRms in Figure IHh- We 
supplement our sample with a compilation of local disks 
(Leroy et al. 2008), as well as with a sample of local 
ULIRGs (Solomon et al. 1997, Rodriguez Zaurfn et al., 
2010). 

Although the actual dependence of SFE on 
sSFR/sSFRms within the MS and the transition 
from MS to the starburst regime will be discussed in 
detail in the next section, here we can derive some 
crucial results about the SFE of the two populations 
(MS and starbursts galaxies) by employing a simple, 
empirical relation between SFE and sSFR/sSFRms- 
Indeed, a linear regression fit and a Spearman's test 
to our high— z data (including starbursts), yields a 
strong and statistically significant correlation (p = 0.71, 
p-value = 0.0012) between star formation efficiency and 
sSFR/sSFRms, with: 

SFE = (8.1 ± 1.2) X [sSFR/sSFRms]^-^'''^"-" (14) 

suggesting substantially higher star formation efficiencies 
for galaxies with enhanced sSFR (Figure |9K). For MS 
galaxies we find an average {SFE) = (14 ± 2) Lo/Mq, 
corresponding to an average gas consumption timescale 
{Tgae = Mgas/SFR) of ~ 0.7 Gyrs, indicative of long- 
lasting star formation activity. This is in direct contrast 
to the short-lived, merger-driven starburst episodes ob- 
served in local ULIRGs, with an average (SFE) — 200 
L0/M0, ~ 12 times higher than that of MS galaxies. 
Since the majority of galaxies at any redshift are MS 
galaxies (e.g., Elbaz et al. 2011) and starbursts galaxies 
seem to play a minor role in the star formation density 
throughout the cosmic time (e.g., Rodrighero et al. 2011, 
Sargent et al. 2012), our results come to support the re- 
cently emerging picture of a secular, long lasting star 
formation as the dominant mode of star formation in the 
history of the Universe. 
To ensure that our result is not an artifact of the as- 
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Fig. 10. — Ratio of direct observables, Lir/Lqq against the 
offset from the main sequence for the same set of objects as in 
Figure 9. The red solid line is the best fit to the data, with a slope 
of 0.80. Coloured tracks corresponds to the scenario by Sargent et 
al. 2012b, following the convention of Figure 9. 

sumed aco values, in Figure [TU] we also plot the ratio of 
the direct observables, Lir/Lqq versus the offset from 
the main sequence for the same set of objects, omitting 
any assumptions for aco- We find the two values to 
strongly correlate with p = 0.79 and a p— value of 3.4 
xl0~^ and a functional relationship of: 

i0.80±0.1 



Lir/L^o = 34.6(±5.0) x [sSFR/sSFRms]' 



(15) 



with MS galaxies exhibiting lower Lir/L^q values by 
a factor of ~ 4. We note, however, that the absence 
of a well defined sample of high— 2; starbursts is strik- 
ing. SMGs, which are frequently regarded as prototypical 
high-z starbursts, have recently turned out to be more 
likely a mixed ensemble of objects, with a significant frac- 
tion of them exhibiting star formation rates and star for- 
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TABLE 3 

Stellar mass, metallicity and gas properties 



Source 


logM* 


.^PP04 " 


log Leo 


«co 




M© 




LqK ^km ^ s pc ^ 


MqK ^km ^ s pc 


ID-8049 


10.97 


8.80 


9.47 


13.4±5.5 


ID-5819 


10.81 


8.76 


9.53 


3.7±1.3 


lD-7o91 


10.83 


8.74 


9.78 


2.4±1.0 


BzK-4171 


10.60 


8.59 


10.32 


2.5±1.1 


BzK-12591 


11.04 


8.63 


10.50 


3.8±2.0 


BzK-25536 


10.51 


8.56 


10.07 


3.1±2.1 


BzK-17999 


10.59 


8.58 


10.23 


3.9±1.2 


BzK- 16000 


10.63 


8.53 


10.32 


9.7±3.1 


BzK-21000 


10.89 


8.60 


10.34 


5.7±2.0 


GN20 


11.36 


9.00 


11.21 


1.0±0.3 


SMM-J2135 


10.34 


8.50 


10.30 


3.2±1.0 


HSLW-01 


10.77 


8.80 


10.30 


1.3±0.4 



Notes: 

a: Derived using the M« — Z relation of Erb et al. (2006), Mannucci et al. 2010 or theoretical arguments (GN20 and HLSW-01). We assume 

a typical uncertainty of ±0.2 for all cases. 



mation efficiencies typical of MS galaxies (Ivison et al. 
2011, Rodighiero et al. 2011). Indeed, while we have only 
three SMGs in this study, we seem to face the same situa- 
tion. Of the three SMGs considered here, only HSLW-01 
appears to be a strong starburst, with sSFR/sSFRms ^ 
10 and SFE - 300 Lq/M©. GN20 is only marginally out- 
side the MS regime, with sSFR/sSFRms ^ 4, although 
with high SFE ~ 120 Lq/Mq, while SMM-J2135 behaves 
like a MS galaxy with sSFR/sSFRms ~ 2.5 and SFE ~ 
30 L0 /M0 , much lower than the star formation efficiency 
observed for local ULIRGs (> 100 Lq/Mq). Clearly a 
larger, objectivcily sc;lc;ctc;d sample of high— 2: starbursts 
is essential to improve this investigation. 

5. VARIATIONS OF SFE WITHIN THE MAIN SEQUENCE 

In the previous section, we demonstrated that the bulk 
of MS galaxies exhibit higher aco values and have lower 
star formation efficiencies than starburst galaxies. Here 
we will attempt to investigate possible variations of SFE 
within the MS, taking into account that the thickness, 
i.e., the spread of the SFR-M* correlation as traced by 
normally star- forming galaxies at any redshift, is not just 
an artifact produced by random noise, but a manifesta- 
tion of the variation of the physical properties of the 
main sequence galaxies, such as color and dumpiness 
(e.g., Salmi et al. 2012; see also Elbaz et al. 2011). Since 
the actual quantity that drives a galaxy above or be- 
low the main sequence is yet unknown, we will consider 
two limiting scenarios, where the relative position of a 
galaxy with respect to the MS is driven by i) variations 
in the gas fraction (/gas = Mgas/[Mgas + M*], or equally 
in the gas to stellar mass ratio, Mgas/M*) while the star 
formation efficiency remains roughly constant within the 
main sequence, or ii) variations in the star formation ef- 
ficiency of MS galaxies, while /gas remains constant, in 
order to explain the observed dispersion in the SFR— M* 
plane. Indeed, in practice there are two ways that a MS 
galaxy could have higher SFR: either it has more raw ma- 
terial (Mgas) to produce stars, or for the same amount 
of Mgas it is more efficient at converting that gas into 
stars. The two scenarios have direct implications on the 
star formation law (SF-law) in the galaxies. The first 
case implies the existence of a global Lm - Mgas relation 
that would apply to all MS galaxies, irrespective of their 
sSFR, like the one presented by Daddi et al. (2010b) and 



Genzel et al. (2010). On the other hand, the second sce- 
nario would imply variations of the star formation law, 
with parallel Mgas - -^'IR relations with constant slope, 
but with a normalisation factor that strongly increases 
with offset from the MS. While a combination of the two 
possibilities might also be plausible, we will consider the 
two limiting cases for simplicity. 

5.1. Limiting Cases 

5.1.1. Scenario I: a Global SF-Law for MS Galaxies; /gas as 
a Key Parameter 

In this scenario, the physical parameter that drives the 
sSFR of a main sequence galaxy is the gas fraction, while 
the star formation efficiency remains roughly constant. 
In this case we have an SF law in the form: 



log Mg 



4as xlogSFR + Ci, 



(16) 



where ^gas = 1 would imply Mgas/SFR — constant. How- 
ever, there is observational evidence that ^gas ~ 0.8, i.e., 
slightly lower than unity (e.g., Daddi et al. 2010, Gen- 
zel et al. 2010, Sargent et al. 2012 in prep.), suggest- 
ing a very mild dependence of SFE from SFR (or from 
sSFR/sSFRms, for fixed stellar mass,) with a slope of 
« 0.2. On the other hand, dividing equation 16 by stel- 
lar mass indicates that Mgas/M* would vary as a function 
of sSFR (or equally as a function of SFR at fixed stellar 
mass) with: 



M, 



gas 



M* 



oc 



sSFR 
sSFRms' 



10.8 



(17) 



In Figure 9b we explore the variations of SFE with nor- 
malised sSFR for three representative bins of mass and 
redshift, as predicted by Sargent et al. (2012 in prep.), 
based on the framework described above. The two fac- 
tors that contribute to the evolutionary tracks are (i) the 
observation that at fixed SFR (or Ljr) starbursts display 
more than an order of magnitude higher SFE than that 
of "normal" galaxies (e.g., Daddi ct al. 2010b, Genzel et 
al. 2010, this study), and (ii) the fact that the depen- 
dence of SFR on Mgas is slightly supra-linear, such that 
SFE increases with Lir even when the main-sequence 
galaxies and starbursts arc considered individually (i.e., 
eq. 16). Point (ii) is reflected in the weak increase of 
SFE throughout the main-sequence and starburst regime 
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(sSFR/(sSFR)ms < 2 and sSFR/(sSFR)ms > 4, respec- 
tively), while the rapid rise in SFE in the transition re- 
gion is the manifestation of point (i) above. Note that the 
"jump" in SFE, from the normal to the starburst galaxies 
is not forced by the assumption of a bimodal star forma- 
tion activity but naturally results from the weak increase 
of SFE with sSFR/sSFRms within the MS. Similar tracks 
are presented in Figure [TOl considering Lir /Lqq instead 

of LiR/Mgas. 

5.1.2. Scenario II: a Varying SF-Law for MS Galaxies; SFE 
as a Key Parameter 

This scenario assumes that /gas for MS galaxies re- 
mains constant, and that the physical parameter that 
dictates the position of the galaxy on the SFR — A/» 
plane is the star formation efficiency. Since, /gas = const, 
then expressing SFE as a function of sSFR/sSFRms, i-e., 
at fixed stellar mass, we have: 



logSFE = log(SFR)-|-C2 



(18) 



suggesting a linear dependence of SFE on sSFR/sSFRms, 
in reasonable agreement with the slope derived by a lin- 
ear fit (^SFE — 1-34 ± 0.13, eq. 14). Although a simi- 
lar slope is found when the fit is redone for MS galax- 
ies alone, excluding the starbursts, a Spearman's test 
suggests that the correlation is statistically insignificant 
{p — 0.58, p— value = 0.13), leaving open the question 
of whether the trend is valid within the MS, or if it is 
driven solely by the starbursts. A similar picture emerges 
when, instead of Mgas, we consider the directly observ- 
able quantity Lqq (Figure [T0| . Again, in this case, ex- 
cluding the starburst data results in a statistically in- 
significant correlation. If indeed the star formation ef- 
ficiency of MS galaxies follows such a steep increase as 
a function of sSFR/sSFRms, aside the existence of par- 
allel Mgas - -^IR SF-laws, it would also imply a smooth, 
continuous transition from MS to starbursts galaxies. 

Both scenarios considered here appear to be consistent 
with our data, and we cannot formally distinguish be- 
tween them based on these individual objects. However, 
as we will see in the next section, the shape of the SED, 
both for individually elected sources as well as for popu- 
lations averaged via stacking, can provide more definite 



DISSECTING THE "THICKNESS' 
SEQUENCE 



OF THE MAIN 



So far, we have seen that the high-quality SEDs of- 
fered by Herschel and mm continuum data can be used 
to reveal a strong dependence of the aco and SEE of 
a galaxy on its deviation from the main sequence (eq.9 
and eq.l4), highlighting significant differences between 
high— z MS galaxies and local ULIRGs, even though the 
former may have (U)LIRG-likc luminosities and SFRs. 
The well-sampled SEDs considered in this study can pro- 
vide the tools to further explore the implications of the 
two scenarios presented in the previous section. 

We have argued that the relation Mgas/Mdust k Z^^, 
observed in the local universe, holds at high— z. A direct 
consequence of this assumption is that the ratio between 
the infrared luminosity and the dust mass of a galaxy at 
any redshift would be: 



Md, 



oc 



^IR 



SFE 



M^ 



X Z 



(19) 



with LiR oc SFR for dusty galaxies where the extinction 
is large. Based on this relation, we can make some pre- 
dictions for the different trends between iiR/Mdust and 
and sSFR/sSFRms, for the two scenarios described 
above. Afterwards, we will attempt to constrain the 
observational trend by stacking large samples of mass- 
selected high-z galaxies. 

6.1. Trends in -/^iR/A/dust for SFE « constant 

For the case where /gas varies within the MS (scenario 
I, §5.1.1), we have Eq. [12] still holding. Also, using the 
z ^ 2 stellar mass - metallicity relation of Erb et al. 
(2006) we get: 

logZ = ez xlogM, + C3 (20) 

with ^ 0.15 for a stellar mass range of 10 < 
log(M*/MQ) < 11. Finally, from the SFR-M* relation, 
we have: 

logSFR = ems xlogM, -)-C4 (21) 
with 0.80. Combining the above relations yields: 

SFE 

log — = (Cms " 4aB * Cms " Cz) X log M* + C5. (22) 

Substituting the values quoted above for the various C 
coefficients very nearly cancels out any dependence on 
logM*, giving: 



, LiR , SFE 

log oc log W C5 

^Mdust ^ Z 



(23) 



i.e., a negligible dependence of the dust-mass- weighted 
luminosity Lm/Mdust with M*. Expressing SFE as a 
function of sSFR, i.e., as a function of SFR for fixed 
stellar mass (and therefore metallicity) yields: 



los 



^IR 



Mdu 



log SFE = (1 - Cgas) X log SFR + Ci, 



(24) 

suggesting a mild dependence of ^ir/ A/dust with sSFR, 
with a slope of ^ 0.2. 

6.2. Trends in iiR/Afdust for /gas ~ constant 

In the case where /gas is constant among MS galaxies 
(scenario II, § 5.1.2), where the thickness of the main se- 
quence is the result of a strongly varying SFE (or equally, 
varying SF-law), we have: 



log 



LiR 
Mdust 



log SFE -I- Ce = log SFR + C7 



(25) 



since Mgas is constant at a fixed stellar mass and /gas. 
This suggests a strong dependence of -Lm/Mdust on 
sSFR/sSFRms- Note that Lm/Mdust is not expected 
to vary with A/*, as was also the case for the previous 
scenario (Eq. [23]). 

Summarizing the analytical predictions, the two sce- 
narios result in two distinct behaviors of LiR/Afdust and 
/gas as a function of sSFR/ sSFRms ■ In the case of a weak 
increase of SFE within the MS (SFE « const, or equally, 
of a global SF-law for MS galaxies), we expect a vari- 
ation of /gas (eq. 17) while iiR/Afdust remains roughly 

In detail, Lip oc (SFRtot — SFRyv), with SFRyv becoming 
significant in the case of low extinction, hence typically at very low 
galaxy masses and metallicities or at high redshifts. 
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Fig. 11. — Top: Observed average SEDs of z ~ 2.0 main sequence galaxies in four sSFR bins, derived by stacking and overlaid with the 
best fit Draine & Li (2007) models. Symbols are color-coded based on the four sSFR bins. The inset panel shows the derived iiR / Afjust a-s 
a function of the offset from the main sequence for the z ~ 2 sample (colored squares) as well as for a set of normal galaxies in the local 
universe by da Cunha et al. 2010a (grey dots). For the z ~ 2 sample, the sSFR/sSFRms is derived based on Herschel measurements. The 
red solid line is the best fit to the local data, with a slope of 0.1. The black solid line shows the linear relation with a 0.5 dex offset. The 
black dashed line shows the expected trend for the case of a steep and continuous increase of SFE with sSFR/sSFRms. Bottom: Observed 
average SEDs of 2 ^ 2.0 main sequence galaxies in three stellar mass bins, derived by stacking and overlaid with the best fit Draine &i Li 
(2007) models. Symbols are color-coded based on the three stellar mass bins. The inset panel shows derived Li^/ M^^^^ as a function of 
the stellar mass for the 2 ~ 2 sample (colored circles) and for the same set of 2 ~ normal galaxies as in the top panel. In both panels the 
SEDs of the stacked samples are normalised to = 1 Lq 



constant within the MS (eq. 25). The opposite trends are 
expected for the case where /gas remains constant (or 
equally, of a varying SF-law), i.e., a steep increase of SFE 
within the MS and a linear increase of iiR,/ M^ust with 
sSFR/sSFRms (eq. 26). We note that both scenarios 
predict only weak dependence of Lm/Mdust with M^,. 
Furthermore, the physical meaning of im/Mdust is the 
luminosity emitted per unit of dust mass, or the strength 
of the mean radiation field heating the dust, and could 
serve as a rough proxy of the effective dust tempera- 
ture. Indeed, since, Lir oc crTd^, two sources with the 
same Lm/Mdust should have similar T^. Similarly, for 
a given Mdust, higher Lir indicates higher Td- There- 
fore, the above analysis also suggests a strong variation 
of the SED shape of the galaxies within the MS for the 
scenario where /gas ~ constant, while only very little, if 
any, change in the shape of the SED for the case where 
SFE PS constant. 

6.3. im/Mdust ond Gas Fractions in MS Galaxies 

Using the stacked samples of 2 ^ 2 MS galaxies in three 
stellar mass and four sSFR bins, described in section 2, 
we can investigate what trends, if any, are present be- 
tween Lm/Mdust- and Lm/Mdust- sSFR/sSFRms- 
To derive the far-IR properties of the stacked samples 
we followed the same SED fitting procedure used for 
the individually detected sources. However, this time we 
take into account the redshift distribution of the stacked 
sources, to account for artificial SED broadening in the 
far-IR and smearing of the spectral features in the mid- 
IR regime. Namely, we construct the whole set of DL07 
models at various rcdshifts and create an average SED 
at the median redshift of the sources considered in the 
stack by co-adding the model SEDs at each redshift, 
weighted by the redshift distribution of stacked sources. 
The two methods return far-IR luminosities that are in 
good agreement, although not accounting for the red- 
shift distribution of the stacked sources results in higher 
.^iR,/.^^dust values. This is expected as in this case the 



artificial broadening drives the fit, erroneously, to higher 
7 values, i.e. higher contribution of the FDR compo- 
nent, that result in lower Mdust for a given Lir. We 
also validate the reliability of the inferred dust masses of 
the stacked samples against a possible artificial broaden- 
ing of the SEDs from the dispersion of the shape of the 
SED of the sources included in the stacking. Namely, 
based on the value of (U) derived from the best fit to the 
stacked data, we generated 1000 artificial SEDs adopting 
a scatter of 0.2 dex in (U) and assigned to each of them a 
redshift based on the redshift distribution of the original 
sample. Then we produced an average SED that we fit in 
the same manner as the original data. We find that the 
derived (U) value for the simulated data is in very good 
agreement with the one obtained for the real data. The 
best fit models for the various sSFR and stellar mass bins 
are shown in Figure[TTl along with the stacked photomet- 
ric points while the inferred parameters are presented in 
Table 4. We note that for every sSFR bin, the derived 
LiR implies an average SFR, and subsequently an aver- 
age sSFR, that is very close to the initial corrected for 
extinction, UV-based sSFR estimate. This is also the 
first time that Herschel data provide evidence, that the 
thickness of the main sequence at z = 2 is real and not a 
due to scatter arising from uncertainties in the derivation 
of SFR. 

As a reference point, we also consider the large compi- 
lation of normal galaxies in the nearby universe ((z) w 
0.05) and local ULIRGs by da Cunha et al. (2010a,b). 
The dust mass estimates in these studies are derived 
based on a two-component fit for warm and cold dust, 
known to give results consistent with those from the 
DL07 models (Magrini et al. 2010). Therefore, a direct 
comparison with our sample is meaningful after applying 
a correction factor of ~ 1.7 to the dust masses of the Da 
Gunha et al. (2010) sample, to account for the different 
K value adopted in their study. 

Figure [TT] reveals a remarkable similarity in the aver- 
age SEDs of z ~ 2 MS galaxies for different sSFR bins, 



18 



10.00 



1.00 



0.10 



0.01 



Mg,./M.«[sSFR/sSFR„s]°'' 
t.„«[sSFR/sSFR„J°' 




^ Normal Galaxit 

Normal Galaxit 



'- 


9 




- 


8 




- 


6 


a) 


- 


4 


s" 


- 


2 ^ 


(d 








- 


1 




- 


05 




, 


01 





10.00 



1.00 



0.10 



0.01 



1 



M /M. « [M.]" 



(X 



^ Z'-~'2 Normal Galaxies 
Z'^O Normal Galaxies 



0.1 



1.0 

sSFR/sSFR 



10.0 



10" 



10'° 10' 
M. [Mo] 



10 



0.9 
0.8 

0.6 
0.4 

0.2 
0.1 
0.05 



0.01 

12 



/gas with increasing stellar mass, with: 

(5.63 ± 0.39) - (0.51 ± 0.10) x log(M, 



Fig. 12. — a) Gas to stellar mass ratios (and gas fractions) as a function of deviation from the main sequence for the stacked sample 
of z ~ 2.0 (filled circles) and the mean values of 2: ~ normal galaxies. Both samples are consistent with Mgas/M* oc [sSFR/sSFRms]"'^ 
suggesting a strong increase of Mg^s/M^ (and /gas) with increasing sSFR at a given redshift. The solid black line is the best fit to the 
z ~ 2 data, and the dashed line has the same slope, normalised to the local sample. The grey horizontal lines define the region of MS 
galaxies, b) Gas to stellar mass ratios (and gas fractions) as a function of stellar mass for the same sample as in the left panel. The data 
suggest a decrease of A/gas/A/* with increasing M, with A/gas/A/, oc A/^'^'^. 

indicating a very small variation of the SED shape of the 
galaxies within the MS. This is further manifested by the 
weak dependence of Lm/Mdust with sSFR/sSFRms (in- 
set panel), both for the z ~ 2 sample and also for the 
local, normal galaxies (slope of 0.11 ±0.04). We reach 
similar conclusions when considering the SEDs for dif- 
ferent stellar mass bins: the shape of the SED and the 
value of Lm/Mdust for both z ^ 2 and z ^ MS galax- 
ies appear to be independent of the stellar mass. The 
observed trends between iin/Mdust vs. sSFR/sSFRms 
and Lin/Mdust vs. are in striking agreement with 
those derived by our analytical approach, favouring a 
small variation of SEE within the MS, and therefore i) 
the existence of a global SE-law for MS galaxies, and ii) 
a step-like dependence of SEE with sSER/sSERms- In- 
deed, the scenario of continuously increasing SEE would 
imply a considerable increase of ii^/Mdust (shown in the 
inset panel of Eigure [TT] (top) with a dashed line) , and 
a noticeable change of SED shape of galaxies within the 
MS, something that is not supported by our data. 

Our data can also be used to infer the dependence 
of /gas on sSFR and M*. Namely, from the derived 
A^dust and Z values, we can use the Sqdr—Z relation 
to estimate A/gas, and subsequently /gas or Mgas/M,. In 
Eigure [T2I we plot the derived Mgas/M, as a function of 
sSER/sSERms and for the normal galaxies at z ^ 2.0 
and z ~ 0. We find a clear trend of increasing Mgas/A^* 
with increasing sSER for both samples, with: 



M, 



M, 



(2.05 ±0.32) X 



sSER 
sSERms' 



i0.87±0.15 



(26) 



very close to the slope derived in equation 17. This re- 
sult further supports that it is variations of /gas 7 rather 
than variations of star formation efficiency, that are re- 
sponsible for the thickness of the SER-A/* relation at 
any redshift. We also reveal a clear trend of decreasing 



(27) 

(M* in Mq) in agreement with various studies that pre- 
dict a similar behavior (e.g., Daddi et al. 2010a, Sain- 
tonge et al. 2011, Dave et al. 2012, Popping et al. 2012, 
Eu et al. 2012). Note that these trends are due to the 
sublinear slopes of the SER vs. and Mgas vs. SER 
relations (see Daddi et al 2010a). 

We recall that the Mgas derived with the (Sgdr method 
traces the total gas of a galaxy, i.e., M^^ + Mm. As 
we have discussed above, in our analysis we have as- 
sumed that H2 dominates the gas mass of relatively 
massive (M* > 10^^ Mq) high-z galaxies. An indi- 
rect way to investigate this assumption is to compare 
the inferred Mgas and Lir of the z ~ 2 sample to the 
L1R-MH2 star formation law for MS galaxies as de- 
rived by Daddi et al. (2010), based on the dynamical 
properties of the individually detected BzK galaxies con- 
sidered in our study. If Mm contributed significantly 
to the total gas mass, then the galaxies would appear 
to show an excess in M^^ with respect to Lir. This 
comparison, for various subsets of our z ^ 2 sample, 
is shown in Eigure 13. The very good agreement be- 
tween the estimated Afgas and the SE-law seems to val- 
idate our assumption (i.e., negligible contribution from 
HI), at least for the stellar mass range considered in this 
study. Indeed, our Mdust measurements are (obviously) 
luminosity- weighted, and thus effectively SER-weighted. 
Similarly, the metallicity estimates that we adopt from 
the literature are luminosity/SER-weighted as well, so 
it is reasonable that our final estimates are most sensi- 
tive to the gas connected to the ongoing SERs, i.e. the 
molecular gas. 

7. THE EVOLUTION OF Lm/Afdust IN MS GALAXIES 

Aside from the similar trend of a weak increase of 
Lm/ Mdust with sSER for z ^ 2 and z ~ normal galax- 
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Fig. 13. — Z/iR vs. Mgas for various stacked subsets in sSFR 
and M« bins (black circles) as well as of the total sample (black 
circle with red dot) of z ~ 2 normal galaxies. The line depicts 
the Lipj-A/h., star formation law of local and high— 2 disk galaxies 
presented by Daddi et al. (2010b). 

ies, Figure [TT] implies that the shape of the SED of MS 
galaxies at a given redshift, as traced by iin/Mdust, is 
not expected to vary significantly with increasing Lir, 
(i.e., SFR). Moreover, normal galaxies at z ~ 2 appear 
to exhibit systematically higher values of LiR/Afdust than 
their local counterparts, suggesting a possible trend of 
LiR/Mdust with cosmic time. To explore this further, we 
enhance our data set by performing SED fitting to the 
mean SED of MS galaxies at z '-^ 1 derived through the 
stacking procedure described in Section 2. The stacked 
data points and best fit models for the z ^ 1 and z ~ 2 
MS galaxies are shown in Figure [TH The final sample 
available for our analysis consists of z ~ normal galax- 
ies from da Cunha et al. 2010 and Draine et al. (2007), 
local ULIRGs from da Cunha et al. 2010b, individually 
detected MS galaxies at z ~ 0.5 and z ~ 1.5, mean SEDs 
of MS galaxies at z ^ 1 and z ^ 2, as well as SMGs at 
various redshifts. We also recall that Lm/Mdust is pro- 
portional to the mean radiation field (U) in the DL07 
models (equation 4). Using the derived (U), Ljr and 
Mdust values for our individuals sources as well as for 
the stacked samples we find an average^* scaling factor 
between the two quantities of Pq ~ 125. 

In Figure[T5](lcft) wc plot Afdust as a function of Ljr for 
the whole set of sources described above. Focusing at 
first on the local samples, we see that z 0, normal 
galaxies follow at tight relation with a slope close to 
but not equal to unity, (^mi = 0.95), possibly mirror- 
ing the iiR — Td relation observed both in the local as 
well as in the high— z universe (e.g., Hwang et al. 2010). 
However, we notice that this is a very mild increase of 
LiR/Mdust- Indeed, the im/Mdust of the whole sample 
of the local normal galaxies exhibits a Gaussian distribu- 
tion with a mean value of (log(LiR/Mdust)} = 2.89 and 
a dispersion of 0.15 dex. To further examine whether 
sources with higher Xir tend to have significantly higher 
-^ir/ A/dust, we split the sample into two luminosity bins. 

There is a small scatter (< 5%) in the derived average 
scaling factor due our definition of LiRas Lg— lOOOpm instead of 
Lm=LQ-oo- 
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Fig. 14. — Observed average SEDs of main sequence galaxies at 
z ~ 1.0 (top) and z ~ 2.0 (bottom) derived by stacking and over- 
laid with the best fit Draine & Li 2007 models taking into account 
the redshift distribution of the stacked sources (grey solid line) and 
assuming the median redshift of each sample (solid black line). The 
black dashed line is the DL07 model without the stellar component 
that is depicted with a green dotted-dashed line. Orange dashed 
and cyan dot-dashed lines show separate contributions of starlight 
and emission from dust heated hy U = Umin (diffuse ISM compo- 
nent) and dust heated by Umin < U < Umax ("PDR" component) 
respectively. The fitted parameters from the best fit Draine & Li 
(2007) model fits are listed within each panel. For comparison the 
best fit SED of the z ~ 1.0, is also shown in the bottom panel with 
a blue-dotted line. 

LiR > 10^^ Lq and Ljr < 10^^ Lq, and find that the two 
sub-samples follow almost identical distributions (Fig- 
ure [T5] left, inset panel). In contrast, local ULIRGs are 
strong outliers from the Afdust-^iR relation formed by 
normal galaxies, with significantly higher Lm^/ M^ust val- 
ues, indicative of their warmer ISM. The IRAS selection 
for the normal galaxies considered here could, in princi- 
ple, introduce a bias towards warmer sources, and miss 
cooler sources whose SEDs peak at longer wavelengths. 
Recent Herschel observations of local galaxies seem to 
confirm this, revealing a population of normal galaxies 
in the local universe that are systematically colder than 
the IRAS selected sample (e.g.. Smith et al. 2012). How- 
ever, while the Herschel-selected galaxies would increase 
the scatter and shift the mean of ^iR/A/dust distribution 
towards a lower value, they are also found to follow a 
LiR — Afdust relation with a slope close to unity, leav- 
ing the overall picture described above unaffected. At 
z = the ULIRGs have an order of magnitude larger 
-^inZ-^^dust than normal (spiral) galaxies, consistent with 
the fact that have a much higher SEE. 

Moving to the high— z samples of MS galaxies, we no- 
tice a small deviation from the local relation, exhibit- 
ing on average higher Lm^/ M^nst values. Having al- 
ready seen that we do not expect large variations of 
Lik/ Mdust within MS galaxies at the same redshift, this 
indicates a small change in the shape of the SEDs of 
MS galaxies towards stronger radiation fields as we move 
back in cosmic time. This is clearly depicted in Fig- 
ure [15] (right) where wc plot Lm/Afdust, or equally (U), 
as a function of redshift, for our samples, along with 
that of local galaxies. Our data suggest an increase of 
(U) oc (1 -I- z)-^'^^ for MS galaxies, pointing also towards 
an increase in their Td with look-back time. This is not 
to be confused with the general conception that ULIRGs 
tend to get colder at higher redshifts. Indeed, when com- 
paring high— z MS ULIRGs to local galaxies of compa- 
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Fig. 15. — Left: Dust mass as a function of infrared luminosity for: local normal galaxies from the sample of da Cunha et al. (2010) 
and the SINGS sample of Draine et al. 2007 (grey points and circles respectively), local ULlRGs from da Cunha et al. 2010b (grey stars), 
individually detected MS galaxies at 2 ~ 0.5 (black squares) and z ~ 1.5 (black squares v^ith blue circle), stacked MS galaxies at z ~ 1.0 
(empty green circle) and z ~ 2.0 (filled green circle) and high— z SMGs (orange stars). The solid line corresponds to a linear relation 
between M^ust and Lir, consistent with the data for local normal galaxies, although the best fit yields a slope of 0.95 (grey dotted line). 
The dashed lines show the linear relation with a ± 0.3 dex offset. Right: Evolution of the mean radiation field (U) oc iiR/M^ust as 
a function of redshift. Black squares with blue circles depict individually detected normal star-forming galaxies at various redshifts. At 
2 ~ we use the average value as derived based on the samples of Draine et al. (2007), da Cunha et al. (2010) and Dale et al. (2012) 
while at high— z the median values from the individually detected sources at z ~ 0.5 and at z ~ 1.5. Black squares represent the stacking 
results at z ~ 1 and z ~ 2. Orange stars represent high— z SMGs considered in the study, while the grey star denotes the position of the 
local ULlRGs based on data from da Cunha et al. (2010b). The black solid line is the best fit to the MS galaxies, yielding a relation of 
((7) oc (1 + z)^'^'^ . The blue and red solid lines show the expected evolution of SFE/Z with redshift for the case of M* = 1 X lO^" Mq and 
= 5.0 X 10"'^'' M0 respectively, highlighting a eo-evolution between SFE/Z and ((7). 



rable iiR (i.e., (U)LIRGs), the latter have much larger 
-^iR/-^^dust values, indicative of a much stronger starlight 
intensity field and a much warmer ISM, in line with re- 
cent studies that find distant ULlRGs to be colder than 
their local counterparts (e.g., Symeonidis et al. 2009, 
Hwang et al. 2010, Magdis et al. 2010, Muzzin et al. 
2010). We also note that the evolution of Td with red- 
shift is directly demonstrated by the comparison of SEDs 
of the z ^ 2 and z ~ 1 samples, with the latter peaking 
at a longer wavelength. Since our stacks consist of galax- 
ies with photometric redshift too, we also explored the 
possibility that the difference in the derived ([/), arises 
from a systematic bias in the redshift distribution of the 
stacked samples. However, we find that for the two sam- 
ples to peak at the same wavelength, this bias should be 
A z w 0.35, which is unlikely given that the systematics 
in the derived photometric redshifts is A z << 0.1 (e.g., 
Grazian et al. 2006, Santini et al. 2009, Pannella et al. 
in prep). 

As we have seen in equation 19, the dust- weighted 
luminosity, (or mean radiation field) is directly propor- 
tional to SFE/Z. It is easy to show how a similar evolu- 
tion with redshift of SFE/Z can be expected based on a 
simple analytical approach. Two trends would drive the 
evolution of SFE/Z. One is the evolution (decrease) of 
metallicity at given stellar mass (Cz(z)) derived by the 
mass-metallicity relation at various redshifts (Tremonti 
et al. 2004 for z=0, Savaglio et al. 2005 for z = 0.5, Erb 
et al. (2006) for z = 2, Sommariva et al. 2012 for z = 3). 
The other is the evolution of the normalization of the 
sSFRms at fixed stellar mass {Cms{.z))i which scales as 



(1- 



At fixed stellar mass, we have: 



lOf 



SEE 



(z) = C + (1 - Cgas) X Cms(z) + Cz(z) (28) 



and the emerging trend of SEE with redshift for the case 
of Af» = 1 X 10^" Mq and = 5.0 x 10^° Mq are shown 
in Figure [15] (right), highlighting the striking agreement 
between the observed and the predicted trends. We 
stress that the metallicity evolution is mapped up to 
z ^ 3, so the values beyond this redshift are based on ex- 
trapolation. However, at very low metallicities, the SER 
and LiR could deviate from their linear relation, with 
LiR, tracing only a small fraction of the total star forma- 
tion activity as dust attenuation becomes small and the 
directly- visible UV emission becomes significant. In that 
scenario, we would expect some fiattening of the SEE 
evolution beyond z ~ 2. In any case, we conclude that 
the observed increase of (C/) with redshift up to z '--^ 2.0 
is in excellent agreement with the expectations based on 
our analytical approach. We stress that the observed in- 
crease of ITJ) with redshift is in line with both alternative 
scenarios, i.e, a weak and strong dependence of SEE on 
sSFR/sSFRms for galaxies within the main sequence. 

7.1. Template SEDs of MS Galaxies 

We have seen that MS galaxies at a given redshift 
tend to have uniform radiation fields, parametrized by 
LiFi/Mdust, and an overall small variation in shape of 
their SED. We also presented evidence that the mean 
radiation field of MS galaxies increases with redshift as 

(l-fz)i-i^ This indicates that we could build template 
SEDs of MS galaxies at different redshifts. In the DL07 
models, the mean radiation field (U) is an interplay be- 
tween 7, which is the fraction of dust emission arising 
from PDRs, and Umin, which is the radiation field heat- 
ing the diffuse ISM. We stress that the same {U) could 
arise from a different combination of 7 and Umin and it 
does not monotonically define the shape of the SED. Al- 
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Fig. 16. — Template SEDs of MS galaxies at various redshift built based on DL07 models and the evolution of (U) with cosmic time as 
(l+z)l l5. The templates assume no evolution of the 7 parameter in the DL07, i.e, the fractional contribution of the PDRs to the total IR 
emission of a galaxy and a small evolution of the qpAH beyond z ~ 2.0, as indicated from our data. The templates are color-coded based 
on their redshift and normalised to ILq. For starburst galaxies we show the best fit template of GN20 (grey line). The black-dashed line 
shows the MS template SED of Elbaz et al. (2011). The templates can be found at http://georgiosmagdis.pbworks.com^ 



though our sample is not sufficiently large to break this 
degeneracy, we can use the mean (U) values at each red- 
shift to construct template SEDs of MS galaxies, under 
the assumption that 7 does not change significantly with 
time. We adopt a constant 7 = 0.02, similar to that re- 
ported for local galaxies by Dale et al. (2012) and Draine 
et al. (2007). Indeed, if we were to reproduce the same 
(U), with higher 7 values then the shape of the SED 
would become unphysically flat close to the peak. Also, 
the galaxies in our sample have a mean 7 very close to 
that value, so it seems like a valid simplification. Then, 
using the relation between (U) and redshift derived in the 
previous section, we built template SEDs of MS galax- 
ies for a grid of redshifts up to z = 2.5 using the DL07 
models. Beyond z > 2.5, we assume a flattening of the 
evolution of (U), similar to that observed for the sSFR 
(e.g., Gonzalez et al. 2010). For the remaining free pa- 
rameter, qpAH, we adopt the mean value of 3.19% up to 
redshift z ^ 1.5 and qpAH = 2.50% at z > 1.5 as indi- 
cated by our data. The resulting MS SEDs are shown in 
Figure [HI 

We note that the overall philosophy is very similar to 
that of Elbaz et al. (2011), who presented a single tem- 
plate for MS galaxies, with a unique far-IR to mid-IR 
luminosity ratio, IRS, at all redshifts up to z > 2 and 
a single template for star-bursting systems with a sig- 
niflcantly higher value of IRS. Here, we extend this ap- 
proach, but introduce some evolution of the SEDs of MS 
galaxies, following the increase of (U) with redshift. Our 
templates result in IRS values that also mildly increase 



from IRS ~ 4.0 at z = to IRS ~ 6.0 at z > 2, in agree- 
ment with Elbaz et al. (2011; and in preparation) and 
with Reddy et al. (2012). Simflar to Elbaz et al. (2011), 
for starbursts galaxies, we show here a unique template, 
using the best flt SED of GN20. Despite the simple and 
straightforward nature of these templates, we note that 
Bethermin et al. (2012b) find that they manage to re- 
produce recent Herschel source counts (e.g., Berta et al. 
2011; Bethermin et al. 2012), including counts measured 
in redshift slices. We note that to investigate the contri- 
bution of an AGN in our templates, we also repeated the 
SED fitting for the stacked samples, excluding this time 
photometric points that are likely to be "polluted" by an 
AGN activity (i.e., Arest < 60 ^m). The inferred (U) val- 
ues are very close to the values derived when considering 
the whole set of photometric points, in agreement with 
various studies that find a negligible effect of an AGN in 
the far-IR colours of a galaxy (e.g., Hatziminaoglou et 
al. 2010, Mullaney et al. 2012). We conclude that our 
templates should be representative for the bulk of MS 
galaxies, even for those that contain an average AGN 
contribution in the mid-IR. The templates can be found 
at http :/ /georgiosmagdis.pbworks. com^ 

8. DISCUSSION 

It seems that the main physical parameter that drives 
the star formation rate of a normal galaxy is the gas frac- 
tion (or Afgas/M*). Indeed, according to our analysis, at 
fixed stellar mass and redshift, there is a strong depen- 
dence of sSFR on Mgas/M*, suggesting that the varia- 
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tions within the main sequence are due to variations in 
the /gas of a galaxy. Interestingly, we find that within the 
MS, Mgas/M* varies almost linearly with sSFR/sSFRms 
(slope of ~0.9). It is important to recall that for the 
gas fractions inferred in this study we use the total gas 
mass of the galaxies. However, Lir, and subsequently 
SFR, is known to strongly correlate with the HCN lu- 
minosity, which is a measure of dense gas {H2 volume 
densities >10'* cm~^), for spirals, local (U)LIRGs and 
QSOs, with an almost linear relation (Gao & Solomon 
2004; Juneau et al. 2009). The inferred dependence of 
SFR to the total gas fraction, along with the fact the wc 
do not observe substantial changes in the star formation 
efficiency, is indicative of a rather constant fraction of 
dense gas among MS galaxies. Alternatively, it suggests 
that the fluctuations of galaxies above and below MS 
are relatively long lasting, allowing the dense gas frac- 
tion to adjust. Indeed, given that the typical timescale 
for dense gas consumption is 10^ yr (and given that it is 
only a fraction of the total gas), fluctuations have to last 
a similar or longer time. This is also the timescale re- 
quired to allow integrated U— V colors to change (Salmi 
et al 2012). Furthermore, the fact that /gas decreases 
with M*, implies that accretion does not keep up with 
SF and outflows at least at z < 2. 

Another important implication of our analysis is that 
any attempt to constrain the evolution of /gas (or 
Mgas/M*) with cosmic time should be focused on nar- 
row stellar mass ranges or take into account the /gas de- 
pendence on both M* and sSFR. Rescaling the /gas es- 
timates of the total stacked samples of 2; ~ 1 and 2; ~ 2 
main sequence galaxies based on the derived trends be- 
tween Mgas/M, - and Mgas/M, - sSFR/sSFRms, 
in Figure 17 we show the evohition of /gas with redshift 
for galaxies of sSFR/sSFRms = 1 and M* = 5 x lO^o 
M© . The evolution of /gas is in excellent agreement with 
the the observationally motivated tracks of Sargent et al. 
(2012b in prep), that are built assuming Mgas oc SFR°-*^ 
at all redshifts and an increase of sSFR as (1 -l-^)^'^ up to 
2; ~ 2.5 followed by a flattening towards higher redshifts. 
We note that recently, Magdis et al. (2012a) provided 
observational evidence for the plateau in the evolution of 
the /gas at z > 2. 

While it seems that there is no value of /gas, stellar 
mass, or sSFR that could be regarded as representa- 
tive for all MS galaxies at a given redshift. we flnd 
that MS galaxies appear to share very similar IR-SED 
shapes, with a characteristic (U) a LiR/M^ust at each 
cosmic epoch. We find that MS galaxies tend to exhibit 
higher Lm/Mdust values as we move back in time, sug- 
gesting that there is a evolution within the MS as a func- 
tion of time towards warmer T^- Nevertheless, we stress 
that high-z MS galaxies with ULIRGs-like Lir, while 
warmer than local, normal galaxies, are colder than local 
ULIRGs, in agreement with various studies that flnd an 
evolution in the of ULIRGs (e.g., Hwang et al. 2010, 
Muzzin et al. 2010). On the other hand, the uniform 
SED shape of MS galaxies at a given redshift suggests 
that Td does not change substantially with Lir within 
the MS. Consequently, it is possible that the strong evo- 
lution of Td with LiR reported by several studies (e.g., 
Magdis et al. 2010, Hwang et al. 2010) is a result of mix- 
ing MS galaxies at various redshifts, or MS and starburst 
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Fig. 17. — Evolution of Mgas/M, with redshift of galaxies with 
sSFR/sSFRms = 1 and Af, = 5 x 10l° Mq. The black squares 
correspond to Mgas/M, of the total stacked samples of z ~ 1 and 
2 ^ 2 nuiin sequence galaxies. The number next to the stacked 
points corresponds to the total number of objects that are included 
in the stack. The open circle at z ~ 0, corresponds to the mean 
/gas of the Leroy et al. (2008) sample. All measurements have 
been corrected to refer to M* = 5 X 10^" Mq. The black dashed 
line depicts the evolution of Mgas / M, with redshift for the case of 
M, = 5.0 X 10^" Mq, as implied by the cosmic evolution of sSFR 
(Sargent et al. 2012 in preparation). 

systems, the latter known to dominate the star-forming 
galaxy population at high luminosities (e.g., Rodighiero 
et al. 2011, Sargent et al. 2012). Although the limited 
number of star-bursting galaxies in this study does not 
allow us to reach solid conclusions, based on Figure 15 
(right) we can speculate that a similar, but less strong, 
evolution of {U), and therefore Td, with look-back time 
is present for this population too (see Bethermin et al. 
2012b for a possible implementation). A larger sample 
galaxies with elevated sSFR with respect to the MS is 
necessary to extend this investigation. 

We stress that the methodology presented in this study 
opens an alternative window for a systematic investiga- 
tion of the gas properties for large samples of high— 2; 
galaxies, for which Herschel has provided a robust char- 
acterisation of their SEDs. While this approach does not 
require CO emission line measurements and is free from 
the uncertain nature of aco, we stress that both tcc;h- 
niques yield consistent results. However, it is important 
to recall the limitations of this methodology. The key 
assumption of this study is that the Sqdj^ — Z relation 
observed in the local Universe does not evolve substan- 
tially with time, implying that roughly half of the metals 
of high— 2: galaxies as well are locked into dust. Evi- 
dently, if the fraction is lower at higher redshifts, then 
the inferred Mgas would be underestimated. However, as 
we argued before, this is unlikely to be a strong effect, 
the opposite might actually be more platisibly expected. 
We also stress that the Sqdr—Z relation is poorly con- 
strained at the low metallicity end {Z < 8.0), where it 
might even become strongly non-linear. However, since 
the lowest metallicity considered in our study is ~ 8.3 
(Figure 4 left, horizontal line), we do not expect to be af- 
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TABLE 4 

Derived properties for stacked sub-samples of z • 



' 2 MS GALAXIES. 





log l^IR 


log JWJugt 


log *^ * ^ 




■^PP04 


log iWgas 




T ^ 


Mq 


Mq 






M0 


sSFRl 


11.44 ± 0.05 


8.18 ± 0.11 


10.09 


4.5±0.7 


8.44 


10.28 ± 0.19 


sSFR2 


11.49 ± 0.05 


8.26 ± 0.09 


10.04 


4.7±0.7 


8.41 


10.38 ± 0.13 


sSFR3 


11.58 ± 0.08 


8.32 ± 0.09 


10.04 


5.2±0.9 


8.41 


10.44 ± 0.16 


sSFR4 


11.80 ± 0.06 


8.51 ± 0.09 


10.17 


5.8±1.0 


8.43 


10.61 ± 0.13 


massl 


11.29 ± 0.06 


8.08 ± 0.10 


9.84 


4.9 ±0.7 


8.35 


10.34 ± 0.15 


mass2 


11.78 ± 0.06 


8.52 ± 0.08 


10.33 


5.0±0.8 


8.48 


10.64 ± 0.13 


mass3 


12.05 ± 0.05 


8.83 ± 0.08 


10.80 


4.9±0.7 


8.65 


10.80 ± 0.13 


Total 


11.55 ± 0.06 


8.27 ± 0.11 


10.10 


5.1±0.6 


8.40 


10.45 ± 0.11 


Total-zl 


11.21 ± 0.09 


8.14 ± 0.13 


10.44 


4.5±0.6 


8.62 


10.04 ± 0.12 



Notes; 

a; The stacking sub-samples in sSFR and stellar mass bins arc given in increasing order or <sSFR>/sS-F-Rms ''■'^d < ^* >• The last line 

corresponds to the total stacked sample of x ~ 1 MS galaxies, 
b: Defined as Lir/Lr. 

c; Estimated through the 5gdr— ^ relation based on the derived Mjuat aud the assumed Zpp()4. 



fected from these uncertainties. Another possible caveat 
is also the assumed metallicities, for which we rely on 
empirical relations and measurements from the litera- 
ture. For example, since the adopted relation to derive 
metallicities for the 2: ^ 2 sample has been calibrated on 
UV-selected galaxies, one could imagine a possible bias 
towards underestimating the actual metallicity of galax- 
ies in our sample (e.g., Onodcra ct al. 2010). If indeed 
the metallicity is underestimated, we would overestimate 
the true Mgas of the galaxies. However, we do not ex- 
pect this effect to impact our main conclusions, given the 
flat shape of the observed stellar mass— Z relation (e.g., 
Erb ct al. 2006). We also note that while stacking data 
for thousands of sources erases interesting details on in- 
dividual sources, it also serves as a tool to average out 
the dispersion in the M* — Z and ^gdr — Z relations. 
However, this method is sensitive to Mjust measurements 
that could, in principle, be affected by systematic biases, 
up to a factor of ~ 3 if one allows for extreme variations 
in the properties of the dust grains. Also, we stress that 
using the templates presented in this study to infer the 
total iiR of individual galaxies from a single detection in 
Rayleigh- Jeans tail, requires extreme caution due to the 
scatter in ([/) which is not fully constrained at high— z. 
Future, direct observations of individual sources in the 
Rayleigh-Jeans tail will guide this investigation and also 
address the physical origin of the scatter, which could 
be linked to the thickness both of the Af, — Z relation 
inside the main sequence as well as of the star formation 
law. Finally, as indicated in Table 4, while the IRS of 
the z ~ 2 stacked samples does not appear to vary with 
stellar mass (hence SFR), there is evidence that it could 
increase with increasing distance to the MS, in agree- 
ment with Elbaz et al. (2011). Variations in IRS, <U > 
and possibly of other physical parameters suggest that a 
single invariant SED for the whole MS should be treated 
as a good flrst order approximation. 

9. CONCLUSIONS 

We have used a sample of 9 normal galaxies at z ~ 
0.5 and z ^ 1.5 and 3 high— z SMGs, all selected to 
benefit from low-J CO observations and high quality far- 
IR photometry, to estimate their gas content and aco 
conversion factors. In addition, large stacked samples of 
> 4000 mass-selected main sequence galaxies at 0.5 < 
z < 2.5 were used to investigate the evolution of their 
Mdust and Mgas content over the last 10 billions years. 



The main findings and implications are summarized as 
follows: 

• For high— z normal galaxies we derive an average 
Q^co value of 5.5 ± 0.4, similar to that of local 
spirals, and ~ 5 times higher that than of local 
galaxies with comparable LiR (i.e., (U)LIRGs). We 
also find that high— z MS galaxies arc ^12 times 
less efficient at converting molecular gas to stars 
when compared to local starbursts, with SEE ^ 14 
L0/M©, corresponding to an average gas consump- 
tion time scale of 0.6 Gyr. 

• For MS galaxies at any redshift, we flnd a weak 

dependence of the dust mass-weighted luminosity 
(ijR/A'/dust): or equally of the mean radiation field 
(?7), on sSFR/sSFRms and M,, complemented by a 
mild evolution with cosmic time as (l-|-z)^-^^. Also, 
for iiR = ig-iooo/iiin, we find average an scaling 
factor between (f/) and iiR/Mjust of Pq ~ 125. 

• We find clear trends of increasing M^^^j Mast (or 
equally of /gas) with increasing sSFR/sSFRmS: 
with a slope of ^/raci ~ 0.9 and decreasing 
Mgas/ Mast with increasing with a slope of 
^frac2 ~ —0.5. Since at a given redshift /gas varies 
with stellar mass and sSFR. any attempt to explore 
the evolution of /gas with redshift should take into 
account these emerging trends. 

• All the above, taken together, lead to the con- 
clusion that variation in the gas fraction (/gas) is 
the driving physical parameter responsible for the 
spread in the SFR— correlation traced by nor- 
mal galaxies at any redshift, and points towards a 
single, tight Lir— Mgas relation for MS galaxies. 

• We find that at a given redshift, MS galaxies seem 
to have uniform IR SED shapes, as parametrized 
by the mean radiation field (JJ) oc LiR/Mjust of 
the DL07 models. Based on this result and on the 
derived evolution of ([/) with redshift, we build and 
provide a set of rcdshift-dcpcndent template SEDs 
of MS galaxies at various redshifts, and show that 
the evolution of the SED of MS galaxies is primarily 
driven by the cosmic evolution of metallicity, and is 
only marginally linked to the rise of average sSFR 
with redshift. 
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Together with the pilot study of Magdis et al. (2011), 
we have demonstrated for the first time the importance 
of studying the continuum SEDs of galaxies in order to 
place constraints on their gas content, through dust mass 
and metallicity estimates. We believe that this is a pow- 
erful technique: until we are in position to formally esti- 
mate aco for high— z galaxies using the resolving power 
of ALMA, it provides an alternative approach for the 
study of Mgas in galaxies throughout out cosmic time. 
However, from our detailed investigation we can point 
out two important results regarding the derivation of 
dust masses, that similar studies should keep in mind; 

• Dust mass estimates based a single temperature 
modified black body model are, on average, a factor 
of ~ 2 lower compared to those derived using the 
physically motivated and more realistic models of 
Draine & Li (2007). 

• While rest-frame submm data (Arcst > 200 fim) re- 
duce the uncertainties in the inferred Mdust by a 
factor of two, excluding mm continuum data from 
the fit does not have an impact on the derived 
Afdust estimates of MS galaxies. On the other hand, 
it appears that in the absence of mm continuum 
data, Mdust estimates of star-bursting systems can 
be grossly overestimated. 
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